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The scientific publications of the United States National 
Museum include two series, Proceedings of the United States 
National Museum and United States National Museum 
Bulletin. 

In these series are published original articles and mono- 
graphs dealing with the collections and work of the Museum 
and setting forth newly acquired facts in the fields of anthro- 
pology, biology, geology, history, and technology. Copies of 
each publication are distributed to libraries and scientific 
organizations and to specialists and others interested in the 
various subjects. 

The Proceedings, begun in 1878, are intended for the 
publication, in separate form, of shorter papers. These are 
gathered in volumes, octavo in size, with the publication date 
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In the Bulletin series, the first of which was issued in 1875, 
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entation. Since 1902, papers relating to the botanical collec- 
tions of the Museum have been published in the Bulletin 
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Two FORTRAN I Programs for the Univariate and Bivariate Analysis 
of Morphometric Data 


Abstract 


Two computer programs, written in FORTRAN II for the 
IBM 7094, perform univariate and bivariate analyses of 
morphometric data and plot bivariate scatter diagrams of 
pairs of data. Basically, the programs are designed to process 
arrays in which rows (up to 100) represent specimens and 
columns (up to 52) represent measured and computed 
variables. Alternatively, an input array may represent a single 
organism, for example a coiled snail or a segmented arthro- 
pod, with the rows representing, respectively, individual 
whorls or segments. A high degree of adaptability to differ- 
ent kinds of problems is achieved because of the numerous 
control cards which specify input and output format, table 
headings, and column numbers of variables on which opera- 
tions are to be performed. 


Program devices allow for the elimination of zeroes rep- 
resenting missing data and for the conversion of selected 
columns to logarithms to the base 10. The computed uni- 
variate statistics include the maximum, minimum, observed 
range, arithmetic mean, standard deviation, standard error 
of the mean, and coefficient of variation. Bivariate statistics 
include the correlation coefficient, slopes and intercepts of 
either regression lines or the reduced major axis, coefficients 
of relative dispersion, and factors for the computation of 
confidence intervals. Finally, the programs plot selected 
pairs of data on bivariate scatter diagrams and furnish the 
end points for the plotting of regression lines or reduced 
major axes. 


Introduction 


Morphometry, or the measurement of form, is becoming in- 
creasingly important in the study of fossil and living orga- 
nisms. The uses and limitations of univariate and bivariate 
techniques in the analysis of morphometric data in paleon- 
tology have been summarized by Imbrie (1956), and the ap- 
plication of such analytic techniques in zoology has been de- 
scribed by Simpson, Roe, and Lewontin (1960). 

In paleontology, it is often necessary to deal with fragmen- 
tary specimens in order to obtain samples (groups of speci- 
mens) which are large enough to be statistically treated. If 
the missing data are represented by zeroes, these zeroes must 
be eliminated from any statistical analysis. Furthermore, if 
the data are transformed to logarithms, a zero representing 
the logarithm of one to the base ten must be treated statis- 
tically and must therefore be differentiated from a zero 
representing missing data. 

Examples of other operations which may be required in 
morphometric work are as follows: 


1. Subtract the thickness of a measuring platform from 
a measurement where, in order to obtain the measurement, 
this thickness had to be included. 

2. Convert certain measurements, which are invari- 
ably in ocular micrometer units, to millimeters. 

3. Convert certain measurements, which are in ocu- 
lar micrometer units only when the specimen is small, to 
millimeters. 

4. Compute the cross-sectional area of a muscle or liga- 
ment, where such an area can be approximated by a sim- 
ple formula, e.g., the area of a circle or a triangle. 

5. Compute measurements which are the sums or dif- 
ferences of other measurements. 

6. Convert any or all measurements to logarithms. 

7. Compute ratios of one measurement to another. 

8. Arrange data in order of decreasing size of one of the 
measurements and/or preserve the order of specimens in the 
input. 
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9. Construct bivariate scatter diagrams of ratios and perform all of these operations. In addition, they punch out 
compute regression lines or lines of best fit. the final measurements and ratios and print out the statistics 
The programs described here, for use on the IBM 7094, listed in Table 1. 


Taste 1.—Symbols and Formulae for Statistics Computed by Data Assembly and Analysis (DASAN) 




















(x or ») = any variate; *see section on statistics) 
Statistics Output FORTRAN = Math Symbol Formula 
Symbol Symbol 
Univariate Statistics 
Number of non-zero values or pairs of values NUM. NUMBER ING, se ee ee ae 
Maximum value MAX. AMAX Dian oY! WE te d= 2 eee eee eee 
Minimum value MIN. AMIN Loca, MO Rea CR fires ae herrea 
Observed range RANGE RANGE OR: OR,=21—22 
vy 
Arithmetic mean MEAN AMEAN z => 
Standard deviation STDDV STDDV = ne V os 
Standard error of the mean SDVMN SDVMN co. a 
=z =z 
100 s 
Coefficient of variation Vv Vv V V= a = 
Bivariate Statistics (Routine #1) = ee 
Correlation coefficient R R r 7 
orrelation coefficien => 
VE(@e—a)'2y—9)* 
Slope of reduced major axis A A a a=2 
7 
Standard error of A SIG(A) STERRA oa o,—4a N 
y-intercept of reduced major axis B B b b=y—awa 
; /2(1—r) (s,?+ sy?) 
Coefficient of relative dispersion from reduced major axis D(RMA) DRMA D D= eee. 100 
Arithmetic mean of » MEAN Y AMEANY 1) en ee ee 
Arithmetic mean of x MEAN X AMEANX ZT) “Al eee Sees S eee eee 
Ordinate of point at which reduced major axis crosses YMAX 1 AVIA 0 en eee ne 
maximum value of x on scale of a plot printed by 
VPLOT 
Ordinate of YMAX 1 in scale-units of a plot printed by TOP 1 TOPIS) ite 1 0) 0) Seater ees e eee eee 
VPLOT 
Ordinate of y-intercept in scale-units of a plot printed BASE 1 TBS Ea Sse ebayer este eee 
by VPLOT 
Bivariate Statistics (Routine #2) 
Correlation coefficient R R r Same as above. 
Z(x— @) (y— 
Slope of regression of » on x A(YX) BYX bie = ne v) 
Y-intercept of regression of » on x B(Y) AY a, ay=Y— bye 
Z(e— 2) (y— 
Slope of regression of x on » A(XY) BXY bs; bay Go v) 
X-intercept of regression of x on B(X) AX az a,=2— bry 
Standard error of estimate in » direction S(Y.X) STESTY Sys gael =m 
Standard error or estimate in x direction SEY) STESTX eee Sy S2V l= 
: : ; 2 100 s 
Coefficient of relative dispersion from regression of » on x D(Y.X) DYX Dre Dyz= 7 
Dyz 





Standard error of D(Y.X) S(DYX) ERRDYX op "pi. V2N 
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Taste 1.—Symbols and Formulae for Statistics Computed by Data Assembly and Analysis (DASAN) (Continued) 
(x or y = any variate; *see section on statistics) 


Statistics Output 
Symbol 
Coefficient of relative dispersion from regression of x on y D(X.Y) 
Standard error of D(X.Y) S(DXY) 
Factor for computing a confidence interval of A(YX)* CONAYX 
Factor for computing a confidence interval of B(Y)* CONBY 
Factor for computing a confidence interval of A(X Y)* CONAXY 
Factor for computing a confidence interval of B(X)* CONBX 


Basically, the programs are designed to process arrays of 
data in which rows represent specimens and columns repre- 
sent measured and computed variables. Alternatively, an in- 
put array may represent a single organism, for example, a 
coiled snail or a segmented anthropod, with the rows repre- 
senting, respectively, individual whorls or segments. A high 
degree of adaptability to different kinds of problems is 
achieved because of the numerous control cards which 
specify input and output format, table headings, and column 
numbers of variables on which operations are to be 
performed. 

The first program, referred to here as DASAN (Data 
Assembly and Analysis), computes the output arrays and all 
statistics; the second program, referred to as VPLOT 
(Variable Plotting), plots the bivariate scatter diagrams, 
using the punched output deck from DASAN as input. 

Background and Acknowledgments.—The initial versions 
of the programs described here were lengthy and were 
designed for the author’s own research concerned with the 
morphometric analysis and phylogenetic interpretation of a 
closely-knit group of species within the Pectinidae (Mol- 
lusca: Bivalvia). Interest expressed by others made clear 








FORTRAN = Math Symbol Formula 
Symbol 
DXY Dey pees 
wr 
D 
ERRDXY a 
“Es °Py 2N 
CONBYX CA, CARS Saas 
Suy N—1 
CONAY CB, CB, = 5. 
VN 
CONBXY CA, CA,= 52 
8. VN— 1 
CONAX CB, CB,= 524. 
VN 


the value of revising the programs so that they would be 
adaptable to a wide variety of problems. 

The reprogramming was carried out at the Smithsonian 
Institution, Washington, D.C., as part of a continuing proj- 
ect concerned with the paleobiology of the Pectinidae, which 
is supported by a grant from the Smithsonian Research 
Foundation. Use of the IBM 7094 of the Columbia Univer- 
sity Computer Center was made possible through the gener- 
osity of Dr. John Imbrie of the Department of Geology. 
Invaluable assistance in constructing the initial programs, 
particularly in incorporating the variable plotting routine 
(Subroutine APLOT), was received from D. M. Vincent 
Manson of the American Museum of Natural History. The 
APLOT subroutine originally appeared in a program written 
by Clarence Bradford and Arthur Gasche of the University 
of Chicago and has subsequently been adapted to a number 
of programs by Manson. Constructive suggestions regarding 
both programs and texts have been gratefully received from 
Alan Cheetham of the Smithsonian Institution, Stephen Jay 
Gould of Harvard University, and Niles Eldredge of 
Columbia University. 


Statistics 


The statistics appearing in the output of DASAN are com- 
puted by means of the formulae shown in Table 1, all of 
which have been taken from either Imbrie (1956) or 
Simpson, Roe, and Lewontin (1960). In addition to their 
use in summarizing data and in the drawing of regression 
lines or lines of best fit directly onto the output of VPLOT, 
the statistics are of use in the construction of modified Dice- 
Leraas diagrams (Simpson et al., p. 355) and in determining 
whether differences in the position and slope of reduced 
major axes are statistically significant (Imbrie, p. 237). 


Using the factors marked by asterisks in Table 1, confi- 
dence intervals (CI) for the slopes and intercepts of regres- 
sion lines may be calculated as follows, where t is the familiar 
Student’s t with N—2 degrees of freedom and other symbols 
are the mathematical symbols of Table 1: 

Cl, =Oyet CAyet 
vz 
CI, =a,+CB,t 
v 
CI, =bey+CArnt 
zy 
Clq =62+CBt 
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DASAN: Description of Program 


Construction and Handling of Arrays.—The construction 
and handling of arrays by DASAN is illustrated in Figure 1. 
Note that the input, referred to as the univariate input array, 
contributes to and becomes part of a univariate output array. 
The latter consists of the univariate input array plus any 
new variables, other than ratios, which are computed from 
the input variables. All of the variables represented in the 
univariate output array may therefore not be univariate 
in a mathematical sense, but they serve here as the data for 
a statistical analysis which yields the univariate statistics 
listed in the introduction. The bivariate array consists only 
of form ratios, which are computed from pairs of variables 
drawn from the univariate output array. 

In each of the three arrays, the first variable (first column) 
consists of specimen identification numbers, with the re- 
maining variables consecutively numbered from left to right 
beginning with variable No. 2. 

In the schematic problem illustrated in Figure 1, the uni- 
variate input array, containing a column of specimen num- 
bers (a) and five columns of measurements (6 through f), 
becomes a subarray (A) within the univariate output array. 
Two of the original variables, d and f, are modified so that 


d 


a 


~ 10.18 
and 
f’=f—9.45 


where 10.18 and 9.45 are constants specified on control cards. 

Four new subarrays have been generated and included 
in the univariate output array. The subarray B contains two 
variables, g and h, where 





= 
2 
and 


h=d/-.,€ 


These are special computations made possible by the inser- 
tion of special statements in DASAN, as described in a fol- 
lowing section. Subarray C consists of variables generated 
by subtracting one of the variables in the preceding columns 
from another. Let us suppose that here 


i=b—e 
and 
j=b—-f 


Subarray D consists of variables generated by adding any 
two of the variables in the preceding columns. For example, 


k=e+h 


Subarray E consists of the logarithms (base 10) of variables 
in the preceding columns. For example, 


1 logiob 


and 
m=logwe 


Univariate statistics for all of the variables (except for 
the specimen identification numbers) in the univariate out- 
put array are computed. Both univariate and bivariate 
statistics are computed for the variables in the bivariate 
array. 

In addition to the bivariate statistics which are computed 
for each variable (ratio) in the bivariate array, an option 
exists whereby a bivariate analysis of x and y transformed 
to logarithms (base 10) may be performed for certain 
columns, the column numbers of which are specified on a 
control card. Let us suppose that in the sample problem 
illustrated in Figure 1 a control card specifies that the 
variables comprising the ratios appearing in columns 2, 3, 
5, 6, and 8 in the bivariate array are to be transformed to 
logs for an additional bivariate analysis. These statistics ap- 
pear in Array IV. The wavy lines beneath the ratios f’/c, 
j/d’, and I/m mark those variables which were not desig- 
nated for logarithmic transformation and for which the 
bivariate statistics from Array III are merely repeated (with 
the exception of the coordinates of the end points of the 
lines of best fit or regressions, which have been altered and 
are no longer valid). 

In the bivariate array for this same sample problem, the 
ratio //m is the ratio of one log (log b) to another (log c). 
Such ratios which are constructed from variables previously 
transformed to logarithms in the univariate output array 
cannot be processed to yield correct univariate and bivariate 
statistics. These invalid statistics are indicated by a wavy line 
in Arrays II and III. Note that correct bivariate statistics 
for this same ratio (log b/log c) appear under b/c in 
Array IV. 

If the horizontal dimension of an output array exceeds 
the width of standard output paper, the array must be 
broken into segments, each of which can be printed on the 
computer paper. This has been done in the sample output 
(see listing of sample output). 

Flow of Control.—The flow of control through DASAN is 
schematically represented in Figure 2, and the main steps 
followed by the program during an execution involving all 
subroutines are listed below: 


DASAN (main) 


1. Control cards, variable formats, and data are read in. 
2. Subroutine RAWTAB is called if so indicated by 


option. 
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Variables—> 
(Measurements) 
albic*'d ie “fF 


ef 
ee 
UNIVARIATE 3¢° 
INPUT ee 5 
ARRAY | e 
7 
8 
; Variables ————> 
aonb. c. dest ghia jy ik Jl 
e}l 
e [2 
UNIVARIATE 2 3 
OUmeUTS ors 
ARRAY e 
7 
8 
Variables (Ratios) —+ 
Ibuich.te'ctetibesbelt 
c dcedd fim 
efi 
Ele 
23 
BIVARIATE a [a 
ARRAY 5 
6 
7 
8 


m 


q 
» 


g 


Variables ———> 
bivchdve ft ig hea yoke ham 


& NUM 
o MAX 
Be eNie UNIVARIATE 
MEAN “ STATISTICS 
STDDV 
SDVMN 
Vv 
Variables —> 
beet: fai js bel 
cdcedd fim 
w” 
2 
3 
5 
wn 





NUM 
MAX 
MIN 
UNIVARIATE 
RANGE 
ean STATISTICS 
STDDV 
SDVMN 
v 


Variables —> 


R 
A 
SIG(A) 

8 BIVARIATE 
(RMA) STATISTICS 
MEAN Y ( y:x) 
MEAN X 

MAX | 
TOP | 
BASE | 


Statistics 








ariables—> 
pic wegiel aoe 
ecdceddfm 

3 R 

o A 

Ss SIG(A) 

a B BIVARIATE 
aoe STATISTICS 
OAM (logy : logx) 

MAX | 
TOP | 
BASE | 


Ficure 1.—Construction of arrays by DASAN for a sample problem. (Abbreviations and construction procedures are explained in the description 


of the program. The bivariate statistics shown here are for the reduced major axis.) 
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Subroutine RAWTAB 


Each of the following steps in RAWTAB is optional, 
depending upon values on control cards. 

3. A read-in constant is subtracted from certain col- 
umns of variables, the column numbers of which are speci- 
fied on control cards. The results replace the original values. 

4. Certain variables, the column numbers of which are 
specified on control cards, are divided by a read-in constant. 
The results replace the original values. 

5. If so indicated by signals on each data card, certain 
data cells, the column numbers of which are specified on 
control cards, are divided by a read-in constant. The results 
replace the original values. 

6. Special variables are computed (e.g., the areas of 
muscle attachments). The results are new columns of vari- 
ables which are added to the original data array. 

7. Certain variables, specified on control cards, are sub- 
tracted from certain other variables, also specified on control 
cards. The results are new columns of variables added to 
the univariate input array. 

8. Same as Step 7, but involving the addition of vari- 
ables. 

9. Certain variables, specified on control cards, are con- 
verted to logarithms to the base 10. The logs are added to 
the univariate input array as columns of new variables. 

10. Control is returned to the main program. 


DASAN (main) 


11. Subroutine XRATIO is called, if so indicated by 
option. 
Subroutine XRATIO 


12. Ratios are computed and entered in a new array. 
The column numbers of variables comprising the ratios are 
specified on control cards. 

13. Control is returned to the main program. 


DASAN (main) 


14. If so indicated by option, subroutine XPUNCH is 
called. 


Subroutine XPUNCH 


15. The univariate output array is punched onto cards 
according to format specified on control cards. 
16. Control is returned to the main program. 


DASAN (main) 


17. If so indicated by option, the first page of the 
univariate output array is printed with no change in the 
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order of specimens in the input deck. 
18. Subroutine UNIVAR is called. 


Subroutine UNIVAR 


19. Univariate statistics are computed from the vari- 
ables which are to appear on one page of output. 

20. The univariate statistics are printed out below the 
page of the univariate output array. 

21. Control is returned to the main program. 


DASAN (main) 


22. If no more univariate output pages remain to be 
printed, the option to proceed to the first page of bivariate 
output is interrogated. If more pages of the univariate 
output array remain to be printed, the next page is printed 
out, followed by a call to UNIVAR and printing out of 
the univariate statistics as in Steps 17 through 21. 

23. If so indicated by option, the first page of the 
bivariate array is printed out, with no change in the order 
of specimens in the input deck. 

24. Subroutine UNIVAR is called. 


Subroutine UNIVAR 


25. Univariate statistics are computed from the bivari- 
ate variables on a page of output. 

26. The univariate statistics are printed out below the 
page of the bivariate array. 

27. Control is returned to the main program. 


DASAN (main) 


28. Subroutine BIVAR is called. 


Subroutine BIVAR 


29. Bivariate statistics are computed for the ratios ap- 
pearing on one page of the bivariate array. 

30. The bivariate statistics are printed out below the 
univariate statistics. 

31. Bivariate statistics for logyy and log,)*. are com- 
puted for any or all of the ratios, as instructed by control 
cards. 

32. The bivariate statistics of the log data are printed 
out below the bivariate statistics of the non-log data. 

33. Control is returned to the main program. 


DASAN (main) 


34. The number of bivariate output pages is checked. 
If more remain to be printed, Steps 23 through 33 are re- 
peated. 
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Read Group I 
Controls 










Read GroupII 
Controls 


Subroutine 
RAWTAB 







Call 
Subroutine 
XRATIO? 









Subroutine 
XRATIO 






Call 
Subroutine 
X PUNCH 2, 

















Punch UV 
and/or BY 


Subroutine 
X PUNCH 





Print univariote 
Statistics 


Print o page of 
UV 


Print bivariate 
log Statistics 


Subroutine 
UNIVAR 












Are 
all UV pages 
printed? 





















Subroutine 
BIVAR 


Convert page 
of BV to logs 









Subroutine 
BIVAR 


Print bivariate 
Statistics 










Are 
all BV pages 
printed 2 






Print univariate 
Statistics 







Order rows of 
UV & BV by 
decreasing x 








UV in order 
of decreasing 







Subroutine 
UNIVAR 






Print one page 
of BY 





data decks to 
process 


No 


Ficure 2.—Simplified block diagram showing the flow of control through DASAN and the major output options. The abbreviations UV and BV 
refer, respectively, to the univariate output array and the bivariate array. 4 and B indicate points between which flow lines cannot be conven- 


iently drawn. 
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35. If no more bivariate output pages remain, the op- 
tion to print all data in order of decreasing values of one 
of the variables, the column number of which is indicated 
on a control card, is interrogated. If so indicated, Steps 17 
through 34 are repeated (flow returns to connecting point 
A in Figure 2) with an option to print the bivariate array 
in order of decreasing values of a specified variable replac- 
ing the option referred to in Step 23. 

36. If more data decks remain to be processed, control 
is passed back to Step 1 (connecting point B in Figure 2) 
with the reading in of new Group-II control cards (see sec- 
tion on input data preparation) and data. Otherwise the 


execution is completed. 


Computation Time.—The test run of DASAN (binary) 
on the IBM 7094, which used the input and furnished the 
output listed and described in following sections, analyzed 
12 specimens, with 28 variables in the input array, 36 vari- 
ables in the univariate output array, and 42 variables in 
the bivariate array, in 52 seconds. 

Core Space.—DASAN in its present form approaches the 
maximum storage allowed in the IBM 7094. Increasing the 
dimensions so that up to 150 rather than only 100 speci- 
mens can be included in a sample resulted in an overlap 
of COMMON and PROGRAM storage areas. The maxi- 
mum capacity of the program (in terms of numbers of speci- 
mens) which can be reached through redimensioning is 
therefore between 100 and 150. 


Listing of FORTRAN II Statements in DASAN 


DASAN (MAIN ROUTINE) 
PROGRAM FOR THE ASSEMBLY AND ANALYSIS OF MORPHOMETRIC DATA 


FEBRUARY 6, 1967 


C 
C 
C 
C THOMAS R. WALLER, DEPT. OF PALEOBIOLOGY, SMITHSONIAN INSTITUTION 
C 
C 
C 


DIMENSION ANAME(21), FMT1(60), FMT2(60) , FMT3(60),FMT4(60),FMT5(60), DMAIN 1 
1FMT6(60) , FMT7(60) , FMT8(96) , IDADD1( 30) , IDADD2( 30) , IDGLAS(30),IDLOG(DMAIN 2 
230) , IDMIN( 30) , IDMTR1( 30) , IEMTR2(5) , IDSUBT( 30), INDXBV(80,4),MMTRNS(DMAIN 3 
3100,5),RATIO(100,56) , RAW(100.53),SPMAX(53),TEMP1(53),TEMP2(56), DMAIN 4 
4 VNAME( 10,24) , SLGMIN( 53) DMAIN 5 

Cc 

COMMON K,RAW,NLOG,NTRNS ,GLASS,OCULAR,NGLASS,IDGLAS,NMMTRN,IDMTR1, DMAIN 6 
LIDMTR2 , NSUB, IDMIN, IDSUBT,NADD, IDADD1,IDADD2,IDLOG,NBIVAR,FMT8,FMT2DMAIN 7 
2, FMT3, FMT6, FMT7, SPMAX, RATIO, INDXBV,MMTRNS , SLGMIN, NBVT DMAIN 8 

Cc 
READ INPUT TAPE 5,34,NDATA,L,NUVTOT,NBVTOT,NRTAB, NRAT DMAIN 9 
READ INPUT TAPE 5,36,NGLASS,NMMTRN,NTRNS ,NSUB, NADD,NLOG DMAIN 10 
IF (NGLASS)30,56,55 DMAIN 11 
55 READ INPUT TAPE 5,36,(IDGLAS(J),J=1,NGLASS) DMAIN 12 
56 IF (NMMTRN)30,58,57 DMAIN 13 
57 READ INPUT TAPE 5,36,(IDMTR1(J),J=1,NMMTRN) DMAIN 14 
58 IF (NTRNS)30,60,59 DMAIN 15 
59 READ INPUT TAPE 5,36,(IDMTR2(J),J=1,NTRNS) DMAIN 16 
60 IF (NSUB)30,62,61 DMAIN 17 
61 READ INPUT TAPE 5,36,(IDMIN(J).J=1,NSUB) DMAIN 18 

READ INPUT TAPE 5,36,(IDSUBT(J),J=1,NSUB) DMAIN 19 
62 IF (NADD)30,64,63 DMAIN 20 
63 READ INPUT TAPE 5,36,(IDADD1(J),J=1,NADD) DMAIN 21 

READ INPUT, TAPE 5, 36,(IDADD2(J),J=1,NADD) DMAIN 22 
64 IF (NLOG)30,66,65 DMAIN 23 
65 READ INPUT TAPE 5,36,(IDLOG(J) ,J=1,NLOG) DMAIN 24 
66 READ INPUT TAPE 5,36,NUVAR,NPAGE1 DMAIN 25 


‘LW 


ORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


READ INPUT TAPE 5,36,NBIVAR,NPAGE2 
IF (NBVTOT)30,68,67 


67 NBVT 
READ 


= NBVTOT - 
INPUT TAPE 


1 
5,36,((INDXBV(I,J),J=1,4),1=1,NBVT) 


68 NPTOT = NPAGE1] + NPAGE2 
READ INPUT TAPE 5,37, ((VNAME(I,J),J=1,24),1=1,NPTOT) 


READ 


1NFMT8 


KK = 
READ 
KK = 
READ 
Kk = 
READ 
KK = 
READ 
KK = 
READ 
KK = 
READ 


Te: 


READ 
READ 
READ 
READ 
READ 
READ 


INPUT TAPE 


12 * NEFMT1 
INPUT TAPE 
12 * NEMT2 
INPUT TAPE 
12 * NFMT3 
INPUT TAPE 
12 * NFMT4 
INPUT TAPE 
12 * NFMT5D 
INPUT TAPE 


12 * NEMT6 
INPUT TAPE 
12 * NFMI7 
INPUT TAPE 
12 * NEMT8 
INPUT TAPE 


INPUT TAPE 
INPUT TAPE 
INPUT TAPE 
INPUT TAPE 
INPUT TAPE 


5, 36,NFMT1,NFMT2 ,NFMT3,NFMT4,NFMT5,NFMT6,NFMT?, 


5,37, (FMT1(1I) ,I=1,KK) 
5,37, (FMT2(I),I=1,KK) 
5,37, (FMT3(1),1=1,KK) 
5,57, (EMT 4(1), 1-1; KK) 


5,37, (FMT5(1I) , I=1,KK) 

5,37 FMI6(I),I - 1,KK) 

537, (FME7(I), I - 1,KK) 

5 37,FMT8(I) ,I=1, KK) 

5, FMT2, (SPMAX(J) J=1,NUVTOT) 

5 yMr2, (SLCMIN(J) ,J-1,NUVTO?) 

5, 30, PUNCH, WRITEL , WRITE? ,RITHT1 , RITHT2 , NORDER 


5 37, ANAME 
5, 39 GLASS , OCULAR, ISIDE,K 


IF (NTRNS)30,46,47 
46 READ INPUT TAPE 5 ,FMT1, ((RAW(I,J),J=1,L),I=1,K 
GO TO 51 
47 READ INPUT TAPE 5 FMEL, ((MMERNS(I,J),J=1,NPRNS), (RAW(I,J),J=1,L), 


1=1,K) 


51 IF (NRTAB)30,49,48 
48 CALL RAWTAB (L,NSPCMP) 


hg Ir (NRAT)30,50,52 
52 CALL XRATIO (NBVT ) 


50 IF (PUNCH) 30,3,2 
© CALL XPUNCH (NUVTOT , NBVTOT , PUNCH) 


3 IF (WRITE1)30,11,4 
4 NPTEMP = 1 
NPTOT = 0 


NV1 = 


2 


Nve = NUVAR 


WRITE] = 0 
5 IF (NPAGEL 


"_ NPTEMP)12,6,6 


10 


6 


ff 
8 
9 


LO 


a 
12 
13 


14 
15 


16 
17 
18 
19 


20 
21 
22 
23 
2h 


a) 
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NPTOT - NPTOT + 1 
WRITE OUTPUT TAPE 6,40, ANAME 


IF (ISIDE)7,9,8 

WRITE OUTPUT TAPE 6,41,K 
GO TO 10 

WRITE OUTPUT TAPE 6,42,K 
GO TO 10 

WRITE OUTPUT TAPE 6,43,K 


WRITE OUTPUT TAPE 6,45, (VNAME(NPTOT,J),J=1,20) 
WRITE OUTPUT TAPE 6,FMT4, (RAW(I,1), (RAW(I,J),J=NV1,NV2),I=1,K) 


CALL UNIVAR (RAW,NV1,NV2,NPTEMP, NUVAR, 1,L,NSPCMP,NSUB,NADD) 
GO TO 5 


NPTOT = NPAGEL 
IF (WRITE2)30,20,13 


NPTEMP = 1 

NVI =e 

NVe = NBIVAR 
NPIMP2 = 1 

NROW1] = 1 

NWROW2 = NBIVAR - 1 
WRITE2 = 0. 


IF (NPAGE2 - NPTEMP)20,15,15 
NPTOr = NPTOT + 1 


WRITE OUTPUT TAPE 6,40, ANAME 


TF ((ISEDE)16, 18 517 

WRITE OUTPUT TAPE 6,41,K 
GO TO 19 

WRITE OUTPUT TAPE 6,42,K 
GO TO 19 

WRITE OUTPUT TAPE 6,43,K 


WRITE OUTPUT TAPE 6,45, (VNAME(NPTOT,J) ,J=1,20) 
WRITE OUTPUT TAPE 6,FMES, (RATIO(I,15, (RATIO(I,J) ,J=NV1,NV2) ,I=1,K) 


CALL UNIVAR (RATIO,NV1,NV2,NPTEMP,NBIVAR,2,L,NSPCMP,NSUB, NADD) 
CALL BIVAR (NROW1 ,NROW2 ,NPIMP2 ) 
GO TO 14 


IF (RITHT1)30,22,21 

WRITEI = 1. 

IF (RITHT2) 30,24 ,23 

WRITE2 = 1. 

ie (RETET + RETHTO)) 30,3225 


Ke =K-1 


DMAIN 73 
DMAIN 74 


DMAIN 75 
DMAIN 76 
DMAIN 77 
DMAIN 78 
DMAIN 79 
DMAIN 80 


DMAIN 81 
DMAIN 82 


DMAIN 83 
DMAIN 84 


DMAIN 85 
DMAIN 86 
DMAIN 87 
DMAIN 88 
DMAIN 89 
DMAIN 90 
DMAIN 91 
DMAIN 92 
DMAIN 93 
DMAIN 9} 
DMAIN 95 


DMAIN 96 


DMAIN 97 
DMAIN 98 
DMAIN 99 
DMAIN1OO 
DMAINLO1 
DMAIN102 


DMAIN1O3 
DMAIN1O4 


DMAIN1O5 
DMAINLO6 
DMAINLOT 


DMAIN108 
DMAIN1O9 
DMATN110 
DMATINL11 
DMAINI12 


DMAIN113 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


DO 29 I=1,Ke2 DMAINI 14 
ee) DMAIN115 
DO 29 Ie=IP1,K DMAIN116 
IF (RAW(I,NORDER) - RAW(I2,NORDER))26,26,29 DMAIN117 
26 DO 27 J=1,NUVIOr DMAIN118 
TEMP1(J) = RAW(I,J) DMAINI19 
RAW(I,J) = RAW(I2,J) DMAIN120 
27 RAW(I2,J) = TEMP1(J) DMAIN121 
DO 28 J=1,NBVTOL DMAIN122 
TEMP2(J) = RATIO(I,J) DMAIN123 
RATTO(L, J) = RATLO(I2, J) DMAIN124 
28 RATIO(I2,J) = TEMP? (J) DMAIN125 
29 CONTINUE DMAIN126 
RECHT = O- DMAIN127 
RITHT2 = 0. DMAIN128 
GO TO 3 DMAIN129 
30 WRITE OUTPUT TAPE 6,44 DMAIN1 30 
GO TO 33 DMAINI 31 
32 NDATA = NDATA - 1 DMAIN1 32 
IF (NDATA)33, 33,70 DMAIN1 33 
70 WRITE OUTPUT TAPE 6,71 DMAIN1 34 
GO TO 1 DMAIN1 35 
33 CALL EXIT DMAIN1 36 
3) FORMAT (I3) DMAIN1 37 
36 FORMAT an DMAIN1 38 
37 FORMAT (12A6) DMAIN1 39 
38 FORMAT (F2.0,4F3.0,13) DMAINL40 
39 FORMAT Re eae e DMAINI41 
4O FORMAT (1H1,21A6) DMAIN142 
41 FORMAT Paseo LEFT VALVES) DMAINL43 
42 FORMAT (1HO,I3,13H RIGHT VALVES) DMAIN1 44 
43 FORMAT rae SPECIMENS ) DMAINI45 
44 FORMAT (1H1,23HANOMALOUS VALVE IN DATA) DMAINIL46 
45 FORMAT vay DMAINL4:7 
71 FORMAT (1H /1H.) DMAINL48 
END 


--------- DASAN (SUBROUTINE RAWTAB) - ----------- 
SUBROUTINE RAWTAB (L,NSPCMP) RAWTB 


DIMENSION DUMMY1( 390) ,DUMMY2(100,56) , IDADD1( 30) , IDADD2( 30) , IDGLAS( RAWTB 
130) , IDLOG( 30) , IDMIN( 30) , IDMTR1( 30) , IDMTR2(5) , IDSUBT( 30) ,MMTRNS( 100RAWTB 
2,5), RAW(100,53) , DUMMY3(80,4) RAWTB 


COMMON K,RAW,NLOG,NTRNS,GLASS ,OCULAR,NGLASS, IDGLAS ,NMMTRN,IDMTR1, RAWTB 
LIDMTR2,NSUB, IDMIN, IDSUBT,NADD, IDADD1, IDADD2, IDLOG,DUMMY1,DUMMY2, RAWTB 


298-840 O—69——_2 


FW 
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2DUMMY 4 , MMTRNS 
DO 28 I=1,K 


IF (NGLASS)30,4,1 
1 DO 3 J=1,NGLASS 
J2 = IDGLAS(J) 
IF (RAW(1I,J2))30,3,2 
RAW(I,J2) = RAW(I,J2) - GLASS 
CONTINUE 


IF (NMMTRN)30,7,5 

DO 6 J=1,NMMTRN 

J2 = IDMTR1(J) 

RAW(I,J2) = RAW(I,J2) / OCULAR 
IF (NTRNS)30,11,8 

DO 10 J=1,NTRNS 

J3 = IDMTR2(J) 

IF (MMTRNS(I,J))9,30,10 

9 RAW(I,J3) = RAW(I,J3) / OCULAR 
10 CONTINUE 


OCOnN] OV On WW Po 


11 J4 =L 
NSPCMP = O 


J4 = J4 + 1 
NSPCMP = NSPCMP + 1 


12 RAW(I,J4) = .5 * RAW(I,21) x RAW(I,20) 


ADD ADDITIONAL COMPUTATIONS LIKE STATEMENT NO. 12 IN PLACE OF 
THESE COMMENTS. EACH COMPUTATION MUST BE PRECEDED BY J4=J4 + 1 AND 
NSPCMP = NSPCMP + 1, AND EACH COMPUTATION STATEMENT BEGINS LIKE 
STATEMENT NO. 12, WITH RAW(I,J4) ON LEFT SIDE OF = SIGN. NO STATE- 
MENT NUMBERS ARE REQUIRED. 


QEATQIATA ATA. QiaAl Ql are 


IF (NSUB)30,18,13 
13 DO 17 J=1,NSUB 

J2 = IDMIN(J) 

J3 = IDSUBT(J) 

einer 

IF (RAW(I,J2))15,14,15 
14 RAW(I,J4) = 0. 

GO TO 17 
15 IF (RAW(I,J3))16,14,16 
16 RAW(I,J4) = RAW(I,J2) - RAW(I,J3) 
17 CONTINUE 


18 IF (NADD)30,24,19 
19 DO 23 J=1,NADD 

J2 = IDADD1(J) 

J% = IDADD2(J) 

J4 = Je +41 


RAWTB 


RAWTB 


RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 


RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 


RAWTB 
RAWTB 


RAWTB 
RAWTB 


RAWTB 


RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWT 

RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 


RAWTB 
RAWTB 
RAWTB 
RAWTB 
RAWTB 


10 
aml 
12 
iy 
14 


i 
16 
a7 
18 
19 
20 
al 
22 
23 
24 


25 
26 


ei 
28 


eo 


30 
mult 
Be 
592 
34 
BD 
36 
a7 
38 
Bo 
4O 


Al 
42 


Aly 
AS 


13 


46 
ei, 
48 
+9 
50 
Sa 


52 
BD 
54 
DD 
56 
D7 
58 
ag 


60 
61 
62 


63 
64 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 
IF (RAW(I,J2))21,20,21 RAWTB 
20 RAW(I,J4) = 0. RAWTB 
GO TO 23 RAWTB 
21 IF (RAW(1I,J3))22,20,22 RAWTB 
22 RAW(I,J4) = RAW(I,J2) + RAW(1I,J3) RAWTB 
23 CONTINUE RAWTB 
24 IF (NLOG)30, 36,25 RAWTB 
25 DO 28 J = 1,NLOG RAWTB 
J2 = IDLOG(J) RAWTB 
Pe J of 1 RAWTB 
IF (RAW(1I,J2))30,28,26 RAWTB 
26 RAW(I,J4) = LOGLOF(RAW(I,J2)) RAWTB 
28 CONTINUE RAWTB 
GO TO 36 RAWTB 
30 WRITE OUTPUT TAPE 6,31 RAWTB 
31 FORMAT (1H1,23HANOMALOUS VALUE IN DATA) RAWTB 
CARI EXIT RAWTB 
36 RETURN RAWTB 
END RAWTB 
NS ah ea) os a DASAN (SUBROUTINE XRATIO) - - ---------- 
SUBROUTINE XRATIO (NBVT) RATIO 
DIMENSION DUMMY1(613), INDXBV( 80,4) ,RATIO(100,56) , RAW(100,53) RATIO 
COMMON K,RAW,DUMMY1, RATIO, INDXBV RATIO 
DO 1 I=1,K RATIO 
RATIO(I,1) = RAW(I,1) RATIO 
Jers RATIO 
DO 1 II=1,NBVT RATIO 
JY = INDXBV(II,1) RATIO 
JX = INDXBV(II,2) RATIO 
Jena sse + 1 RATIO 
1 RATIO(I,J2) = RAW(I,JY) / RAW(I,JX) RATIO 
RETURN RATIO 
END RATIO 
eae oe a DASAN (SUBROUTINE XPUNCH) ------------ 
SUBROUTINE XPUNCH (NUVTOT,NBVTOT, PUNCH) PUNCH 
DIMENSION DUMMY1 (320) ,DUMMY2(173) , FMT2(60) , FMT3(60) ,RATIO(100, 56) , PUNCH 
LRAW(100,53) Bes 
COMMON K,RAW,DUMMY1,FMT2,FMT3,DUMMY2, RATIO PUNCH 
Bp CPUNGH = 2.1.2, 1 PUNCH 


1 PUNCH = PUNCH - 1. PUNCH 


i 
HOW ONAN F WW 


PH 
Wf 
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G 
WRITE OUTPUT TAPE 7,FMT2,((RAW(I,J),J=1,NUVTOT),I=1,K) PUNCH 
C 
TF (PUNCH = 1.)3,3,2 PUNCH 
2 WRITE OUTPUT TAPE 7,FMT3,((RATIO(I,J),J=1 NBVTOT),I=1,K) PUNCH 
c 
3 RETURN PUNCH 
END 
Chowan te Sue ose DASAN (SUBROUTINE UNIVAR) --------+---- 
C 
SUBROUTINE UNIVAR (XX,N1,N2,NPTEMP,NVAR,M,L,NSPCMP, NSUB, NADD) UV 
C 
DIMENSION AMAX(20) , AMEAN( 20) , AMIN(20) , ANUM(20) , BNUM( 20) , DIFFER( 100UV 
1,20) ,DIFSQR(100,20) , DUMMY1(100,53) , DUMMY2(192) , FMT6(60) ,FMT7(60), UV 
2TDLOG2(100,20) , NUMBER(20) , RANGE(20) ,SDVMN(20) ,STDDV(20) ,SUM(20), UV 
3SUMSQR( 20) ,V(20) ,X(100,20) ,XX(100, 56) , IDLOG(30) , DUMMY3(217) UV 
C 
COMMON K,DUMMY1,NLOG,DUMMY2 , IDLOG, DUMMY3, FMT6, FMT? UV 
é 
po 43 I = 1,K UV 
DO 43 J = 1,NVAR UV 
MZ IDLOGe (had), =O UV 
IF (M - 1)28,44,29 UV 
44 TF (NLOG) 34,29, 38 UV 
38 DO 33 Jl=1,NLOG UV 
J2 = L + NSPCMP + NSUB + NADD + Jl UV 
IF (J2-N1)33,30,30 UV 
30 IF (J2-N2)31,31,29 UV 
31 J3 = IDLOG(J1) UV 
cy = J2-NP 4 1 UV 
DO. 32) Tai K UV 
TP CXKCE IS S45 52 5.55 UV 
35 IDLOG2(1,J4) = 1 UV 
32 CONTINUE UV 
33 CONTINUE UV 
C 
29 DO 1 J=N1,N2 UV 
an, af es ay es al UV 
DO 1 I=1,K UV 
PMG od) =o XC Ta) UV 
NN1 = 1 UV 
NN2 = N2 - Nl i+ 1 UV 
C 
DO 3 J=NN1,NN2 UV 
AMAX(J) = X(1,3) uM 
DO 3 I = 2,K UV 
IF (AMAX(J) - X(I,J))2,3,3 UV 
2 AMAX CS) = XCI yd) UV 
3 CONTINUE UV 
C 
DO 14 J=NN1,NN2 UV 
It Sal UV 


19 


hb 


ON UW FW 


10 
11 
12 
13 
14 
ID 
16 
LZ 
18 
1g 
20 
en 
22 


2D 
24 
25 
26 
27 
28 


ag 


31 
32 
33 
34 


35 
36 


WI wo 


NI 


40 
10 
iP 
2. 


13 
14 


1D) 
16 


a 
18 


aug) 


20 


4a 
22 
235 


24 
a2 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


ANUM(J) = O. 

SUMCS) "= O% 

IF (IDLOG2(1I,J))28,4,7 
ROR CEA) |) 76557 

IF (I-K)6,7,28 

peters teen 

GO TO 39 

AMIN(J) = X(I,J) 

SUMCIe) o=— eC Ladi) 

ANUM(J) = 1. 

ies tee a 

Tih (i =.K 40,4013 

IF (IDLOG2(1I,J))28,9,10 
my (XC Tt) )VO.64 VO 
ANUM(J) = ANUM(J) + 1. 
SUM(J) = SUM(J) + X(I,J) 
i CAMENGH) “= X(T) )8, 8,11 
AMIN(J) = X(I,J) 


GO TO 8 
RANGE(J) = AMAX(J) - AMIN(J) 
AMEAN(J) = SUM(J) / ANUM(J) 


DO 19 J=NN1,NN2 
Tee, tll 

SUMSQR(J) = 0. 

GE CUDLOGA(iad ) eo. 15>, 17 

Te ONC T yy) 7G 7. 

a totee 

Tie Chak) LO. 16, 19 

DIFFER(I,J) = X(I,J) - AMEAN(J) 
DIFSQR(I,J) = DIFFER(I,J) ** 2. 
SUMSQR(J) = SUMSQR(J) + DIFSQR(I,J) 
GO TO 16 

CONTINUE 


DO 20 J=NN1,NN2 
STDDV(J) = (SUMSQR(J) / (ANUM(J) - 1.)) ax -5 
V(J) = 100. « STDDV(J) / AMEAN(J) 


DO 23 J=NN1,NN2 
NUMBER(J) = 0 

DOpe2 al 1K 

IF (IDLOG2(1,J))28,42,21 
TCO I) ete. 

NUMBER(J) = NUMBER(J) +1 

CONTINUE 

BNUM(J) = NUMBER(J) 

SDVMN(J) = STDDV(J) / BNUM(J) «« .5 


GO TO (25,26) ,M 


qaaaq 


agaaq 


a9gee4aq geaeggdesaaa 


UV 


UV 


WRITE OUTPUT TAPE 6,FMT6, (NUMBER(J),J=NN1,NN2), (AMAX(J) ,J=NN1, NN2)UV 


37 
38 
39 
ho 
41 
42 
43 
hy 
45 
Lé 
re) 
4g 
4g 
50 


be 
2D 
54 
5D 
56 
Dd 


58 
29 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 


a 
73 


7 
a2 
76 
77 
78 
(és) 
80 
81 
82 


83 
8h 
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1, (AMIN( J) ,J=NN1,NN2), (RANGE(J),J=NN1,NN2) , (AMEAN(J) , J=NN1,NN2) , (STUV 


2DDV(J),J=NN1,NN2), (SDVMN(J),J=NN1,NN2), (V(J),J=NN1, NN2) 
GO TO 27 


UV 
UV 


26 WRITE OUTPUT TAPE 6,FMT7,(NUMBER(J),J=NN1,NN2), (AMAX(J) ,J=NN1,NN2)UV 
1, (AMIN(J),J=NN1,NN2) , (RANGE(J) ,J=NN1,NN2) ,(AMEAN(J) , J=NN1,NN2) , (STUV 


34 


ei, 


36 
28 


4,2) ,X(100,53) , YMAX1(20) , LGCOL(20) , DUMMY4(100,5) ,SLGMIN(53) 


52 
22 


2D 
34 


33 


2DDV(J),J=NN1,NN2), (SDVMN(J),J=NN1,NN2),(V( J), J=NN1, NN2) 
GO TO 27 


WRITE OUTPUT TAPE 6,36 
CALL EXIT 


NPTEMP = NPTEMP + 1 

Nl = N2 + 1 

N2 = Ne + NVAR - 1 

FORMAT (1H1,23HANOMALOUS VALUE IN DATA 


RETURN 
END 


Seay ee DASAN (SUBROUTINES B EVR): 9 eset ee 


SUBROUTINE BIVAR (NROW1,NROW2, NPTMP2 ) 


UV 
UV 


UV 
UV 


UV 
UV 
UV 
UV 


UV 
UV 


BVRMA 


DIMENSION A(20),AMEANX(20) , AMEANY(20) ,B(20), BASE1 (20) , DIFFER(100,2BVRMA 
1) ,DIFSQR(100,2) ,DRMA(20) , DUMMY1 (222) , DUMMY2 (240) , DUMMY3(100,56) , FMBVRMA 
2T8(96) , INDXBV(80,4) , PAIRS(20) ,R(20) ,RSQ(20) ,SPMAX(53) ,STDDVX(20), BVRMA 
3STDDVY(20) ,STERRA(20) ,SUM2(2) , SUMMUL( 20) , SUMSQ2(2) , TOP1 (20) , TX( 1LOOBVRMA 


BVRMA 


COMMON K,X,NLOG,DUMMY1,NBIVAR,FMT8 ,DUMMY2,SPMAX, DUMMY3, INDXBV, DUMMBVRMA 


1Y4,SLGMIN,NBVT 


M ="0 


3 


0 
O 

DO 16 II = NROW1,NROW2 
IF (II - NBVT)55,55,54 


dl = ENDXBVCIT.1) 
J2 = INDXBV(II,2) 
LL = INDXBV(II,3) 


LGT = INDXBV(II,4) 

TCM) 22555. 55 

IF (INDXBV(II,4))34,33, 34 
IM = IM+1 
LGCOL(IM) = LL 
PAIRS(LL) = 0. 


SUMSQ2( 1) = 0. 
SUMSQ@2( 2) = O. 
SUM2( 1) = 0. 
SuUM2( 2) = 0. 
SUMMUL(LL) = 0. 


BVRMA 


BVRMA 
BVRMA 
BVRMA 
BV RMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BV RMA 
BVRMA 
BV RMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 


85 
86 
87 


88 
89 
90 
gl 


92 
93 


94 
93 
96 
97 


98 
99 


OW FW Po = 


oon) 


10 
ane 
2 
13 
14 
15 
16 
1? 
18 
ee) 
20 
2a 
22 
aD 
24 
25 
26 
ey, 


41 


29 


30 


24 
25 


60 
pit 


61 


api: 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


IF (M)30,41,30 
tf (GT) 29.30.29 
LGT = 0 

MMig—srele 


Title 

DO 5 I=1,K 

Tey AGG) 2. 2 

xed) ei. dil) 

TCE 2) = XC J2) 

MeV CKL sot) 20,5, e4 

OK CE 2 2085.25 

MGhad) = LOGLORC XKCD, Jl)!) 
KEE Je) = LOGLOM (X(1,,J2)) 

me GZ) 4.60.4 

IF (SPMAX(J1))20,20,51 
TSPMAX = SPMAX(J1) 

SPMAX(J1) = LOGLOF(SPMAX(J1)) 
IF (SPMAX(J2))20,20,61 
TSPMXX = SPMAX(J2) 

SPMAX(J2) = LOG1OF(SPMAX(J2) ) 
Tie es 

GO TO 4 

MC Chad) ) 5.5.3 

be GOR Cl Je) as >t 

PAIRS(LL) = PAIRS(LL) + 1. 
SUM2( 1) = SUM2( 1) + X(I,J1) 


SUM2( 2) = SUM2( 2) + X(I,J2) 
CONTINUE 

AMEANY(LL) = SUM2( 1) / PATIRS(LL) 
AMEANX(LL) = SUM2( 2) / PAIRS(LI,) 
DO 6 J=1,2 

DO 6 I=1,K 

DIFSQR(I,J) = 0. 

DIFFER(I,J) = 0. 


DOL 9. © = 1,K 
ie (LGT)47,7,47 

rH CTIKCH 1) 20,9 ,25 

EC PK Gl... 2))i205,9';,10 

iE WGx GL dil) 91,0 

iS (Chad) ) 10.9, 10 

DIFFER(I,1) X(I,J1) - AMEANY(LL) 
DIFFER(I,2) X(I,J2) - AMEANX(LL) 
DIFSQR(I,1) DEPEERC 1) ae 12 
DIFSOR(1,2) = DIFFER(L,2) ** 2. 
CONTINUE 


DO V4 T=1 ,K 
DIFMUL = DIFFER(I,1) * DIFFER(I,2) 


BVRMA 
BVRMA 
BVRMA 
BVRMA 


BVRMA 


BVRMA 
BVRMA 
BVRMA 
BV RMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 


BVRMA 
BV RMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 


BVRMA 
BVRMA 


BVRMA 
BVRMA 
BVRMA 
BVRMA 


BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 
BVRMA 


BVRMA 
BVRMA 


17 


28 
ao 
30 
Fil 


32 
33 
3 
35 
36 
37 
38 
39 
ho 
hy 
ke 
43 
bik 
45 
h6 
47 
48 
49 
50 
Silk 
52 
53 
54 
55 


Di 
58 


9 
60 
61 


62 


644 
65 
66 
67 
68 
69 
70 


72 


73 
74 


UNITED STATES NATIONAL MUSEUM BULLETIN 285 








SUMMUL(LL) = SUMMUL(LL) + DIFMUL BVRMA 75 

IF (LGT)26,27,26 BVRMA 76 
265mm (TXCL.1) 20,1428 BVRMA 77 
28 obE (CLKGh.2))20. aa 1s BVRMA 78 
27 TR OC a ke eee BVRMA 79 
2 CRG eesop less dees BVRMA 80 
13 SUMSQ2(1) = SUMSQ2(1) + DIFSQR(I,1) BVRMA 81 
SUMSQ2( 2) = SUMSQ2( 2) + DIFSQR(I, 2) BVRMA 82 

14 CONTINUE BVRMA 83 
STDDVY(LL) = (SUMSQ2(1) / (PAIRS(LL) - 1.)) «+ .5 BVRMA 84 
STDDVX(LL) = (SUMSQ2(2) / (PAIRS(LL) - 1.)) ** .5 BVRMA 85 
SY2 = STDDVY(LL) +#* 2. BVRMA 86 
SX2 = STDDVX(LL) «** 2. BVRMA 87 
R(LL) = SUMMUL(LL) / (SUMSQ2( 1) * SUMSQ2( 2)) «x .5 BVRMA 88 
AC(LL) = STDDVY(LL) / STDDVX(LL) BVRMA 89 
RSQ(LL) = R(LL) «« 2. BVRMA 90 
STERRA(CEL). =) ACh) « CC. = RSQCLE)) / PATRS GEE)! Sn) -5 BVRMA 91 
B(LL) = AMEANY(LL) - AMEANX(LL) »* A(LL) BVRMA 92 
DRMA(LL) = 100. * (2. * (1.-R(LL)) * (SX2+SY2) / ((AMEANX(LL)**2.)BVRMA 93 

1 + (AMEANY(LL)«+*2.))) ** .5 BVRMA 94 
YMAX1(LL) = SPMAX(J2) * A(LL) + B(LL) BVRMA 95 

IF (LGT)49,48,49 BVRMA 96 

48 VDY = SPMAX(J1) / 49. BVRMA 97 
GO TO 50 BVRMA 98 

4Q VDY = (SPMAX(J1) - SLGMIN(J1)) / 49. BVRMA 99 
50 TOP1(LL) = (YMAX1(LL) - SPMAX(J1)) / VDY BVRMA100 
BASE1(LL) = B(LL) / VDY BVRMA101 

DO 17 I=1,K BVRMA102 

TS CLED ase NGn > BVRMA103 

USP KCL, Jit) = TX CE.) BVRMA104 
7 eCh 2) = 2x2) BVRMA105 
SPMAX(J1) = TSPMAX BVRMA106 
SPMAX(J2) = TSPMXX BVRMA107 

16 CONTINUE BVRMA108 
54 NN = NBIVAR BVRMA109 
IF (M+MM) 36,31, 36 BVRMA110 

31 IF (89-K)45,46,46 BVRMA111 
45 WRITE OUTPUT TAPE 6,19 BVRMA112 
46 WRITE OUTPUT TAPE 6,18 BVRMA113 
GO TO 39 BVRMA114 

36 IF (M)43,42,43 BVRMA115 
42 IF (70-K)44,46,46 BVRMA116 
44 WRITE OUTPUT TAPE 6,19 BVRMA117 
43 WRITE OUTPUT TAPE 6,40, (LGCOL(J),J=1, IM) BVRMA118 


39 WRITE OUTPUT TAPE 6,FMT8,(R(J),J=2,NN),(A(J),J=2,NN) , (STERRA(J) , J=BVRMA119 
12,NN),(B(J),J=2,NN) , (DRMA(J),J=2,NN), (AMEANY(J),J=2,NN) , (AMEANX(J) BVRMA120 
2,J=2,NN), (YMAX1(J),J=2,NN), (TOP1(J) , J=2, NN), (BASE1(J) ,J=2, NN) BVRMA121 


DO 53 J=2,NN BVRMA122 
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R(J) Oz 
A(J) O. 
STERRA(J) 
Budo) =O: 
DRMA(J) = 
AMEANY (J) O. 
AMEANX(J) O. 
MMAXACS) =O: 
TOPICS) = 0. 
53 BASE1(J) = 0. 


| 
oO 


eS: 


IF (MM)37,38,37 
37 M=1 
MM = O 
GO TO 32 
38 NPTMP2 = NPTMP2 + 1 
NROW1 = NROW2 + 1 
NROW2 = NROW2 + NBIVAR - 1 
GO TO 22 


20 WRITE OUTPUT TAPE 6,21 
CALL EXIT 


21 FORMAT(1H1,23HANOMALOUS VALUE IN DATA) 


£9 


BVRMA123 
BVRMA124 
BVRMA125 
BVRMA126 
BVRMA127 
BVRMA128 
BVRMA129 
BVRMA130 
BVRMA131 
BVRMA132 


BVRMA133 
BVRMA134 
BVRMA135 
BVRMA136 
BVRMA137 
BVRMA138 
BVRMA139 
BVRMA140 


BVRMA141 
BVRMA142 


BVRMA143 


18 FORMAT (1H-,93HBIVARIATE ANALYSIS (REDUCED MAJOR AXIS) OF Y/X. REBVRMA144 


1SULTS NOT VALID FOR COLUMNS OF LOG Y/LOGX.) 
19 FORMAT (1H1) 


BVRMA145 
BVRMA146 


4O FORMAT (1H-,80HBIVARIATE ANALYSIS (REDUCED MAJOR AXIS) OF LOG Y/LOBVRMA147 


1G X VALID ONLY FOR COLUMN(S) ,10I4/1H,10L4) 


22 RETURN 
END 


Bee ee Soto ea jo DASA (CSUSROULENE: BEVAR)ie =e a0 =) === ete on 
ALTERNATE SUBROUTINE BIVAR WHICH COMPUTES 
THE REGRESSIONS OF Y ON X AND X ONY. 


SUBROUTINE BIVAR (NROW1, NROW2 , NPTMP2 ) 


DIMENSION AMEANX(20),AMEANY(20),AX(20) ,AY(20) , BASE2(20) , BASE3(20), 
1BXY(20) , BYX(20) , CONAX(20) , CONAY(20) , CONBXY( 20) , CONBYX(20) , DIFFER(1 
200,2) ,DIFSQR(100,2) , DUMMY1( 222) , DUMMY2( 240) , DUMMY3(100, 56) , DUMMY4( 
3100,5),DXY(20) , DYX(20) , ERRDXY(20) , ERRDYX(20) , FMT8(96), INDXBV( 80,4 
4) ,LGCOL(20) , PAIRS(20) ,R(20) , RSQ(20) , SLGMIN(53) , SPMAX(53) , STDDVX(20 
5), STDDVY(20) , STESTX(20) , STESTY (20) , SUM2(2) , SUMMUL(20) , SUMSQ2(2) , SX 
6Y(20) ,SYX(20) , TOP2(20) , TOP3(20) , TX(100,2) ,X( 100,53), XMAX(20) , YMAX2 
7(20) 


COMMON K,X,NLOG,DUMMY1, NBIVAR, FMT8 , DUMMY2 , SPMAX , DUMMY3, INDXBV, DUMM 
1¥4,SLGMIN, NBVT 


Mie 


BVRMA148 


BVRMA149 


BVYXY 


he 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


O CON DU FW Pf 


BVYXY 10 
BVYXY 11 


BVYXY 12 


20 


32 
D2 


D5 
34 


33 


41 
29 


30 


24 
29) 


60 
Dal 


be 


FW fo 
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0 
0 

DO 16 II = NROW1, NROW2 
IF (II - NBVT)55,55,54 
Jl = INDXBV(II,1) 

J2 = INDXBV(II,2) 

LL = INDXBV(II,3) 

LGT = INDXBV(II,4) 

in (M22 33435 

IF (INDXBV(II,4))34,33, 34 
IM = IM+1 

LGCOL(IM) = LL 
PAIRS(LL) = 0. 


2 


SuMSQ2( 1) = 0. 
SuMSOg2( 2).= 0. 
SUM2( 1) = 0. 
SUM2( 2) = 0. 
SUMMUL(LL) = 0. 


IF (M)30,41, 30 
IF (LGT)29,30,29 


Gh =) 0 
MM = 1 
Z=0. 
DOORS EK 


GEG) 2. ab 

Glee) = Cg) 

SG. 2) aera 2) 

EE OCXGEZI1)) 20), 5,24 

EF (kCEsd2))20,.5. 23 

MCh dl) = hOG1LOr Cuts s)) 
XC, J2)) = LOGLor (x(1,.J2)) 

IF (Z) 4,60,4 

IF (SPMAX(J1))20,20,51 
TSPMAX = SPMAX(J1) 

SPMAX(J1) = LOG1LOF(SPMAX(J1)) 
IF (SPMAX(J2))20,20,52 

TSPMXX = SPMAX(J2) 

SPMAX(J2) = LOG1OF(SPMAX(J2)) 
Tiga lee 

GO TO 4 

TE Ces dl) 3.15). 3 

DE CXC 2) ye 5. 4 

PATRS (iL) = PATRS (iL). 4.1". 


SUM2(> 1) = sSUM2( 1) + XG J1) 
SUM2C 2) s="SUM2( 2) Fe axCis2) 
CONTINUE 

AMEANY(LL) = SUM2( 1) / PAIRS(LL) 
AMEANX(LL) = SUM2( 2) / PAIRS(LL) 
DO 6 J=1,2 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 


BVYXY 


ILS) 
14 
15) 
16 
aa 
18 
19 
20 
2a 
22 
25 
a4 
25 
26 
Zn 
28 
29 
30 


31 
32 
BD 
34 


35 
36 
37 
38 
39 
LO 
41 
42 
43 
iy 
45 
46 
47 
48 
49 
50 
5 
52 
53 
54 
55 
56 
57 
58 


29 
60 


61 


4? 
25 


10 


alia 


26 
28 
All 
a2 
13 


14 
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DOnowl=1s, Ke 
DIFSQR(I,J) = 0. 
DIFFER(I,J) = O. 


DOw9 y=) LK 

Ir (LGT)47,7,47 
TRACT) 205 9429 

Tr (TMC. 2))20), 9,10 

ie Cx(iedl))8,9,0 

Ti OXGiewde)») 10),.9), LO 

DIFFER(I,1) = X(I,d1) - AMEANY(LL) 
DIFFER(1I,2) = X(I,d2) - AMEANX(LL) 
DIFSQR(I,1) = DIFFER(I,1) ** 2. 
DIFSQR(I,2) = DIFFER(I,2) ** 2. 
CONTINUE 


DO ms: =], K 

DIFMUL = DIFFER(I,1) + DIFFER(I,2) 
SUMMUL(LL) = SUMMUL(LL) + DIFMUL 
IF (LGT)26,27,26 

Tm (TX(L, 1) 20), 14,28 

Tia (mMGie 2) 2Ojd4., 19 

Tr Ou Wada.) alia, a a2 

ir (XCLed2) )i5514),.15 

SUMSQ@2(1) = SUMSQ2(1) + DIFSQR(I,1) 
SUMSQ2(2) = SUMSQ2(2) + DIFSQR(I,2) 


CONTINUE 
STDDVY(LL) = (SUMSQ2(1) / (PATRS (LL) = 2.)) e# <9 
STDDVX(LL) = (SUMSQ2(2) / (PATRS(L)) = 1.) + #8 eo 


BYX(LL) = SUMMUL(LL) / SUMSQ2(2) 
BXY(LL) = SUMMUL(LL) / SUMSQ2(1) 


AY(LL) = AMEANY(LL) - BYX(LL) + AMEANX(LL) 
AX(LL) = AMEANX(LL) - BXY(LL) * AMEANY(LL) 
sy2 STDDVY(LL) ** 2. 


See STDDVX(LL) ** 2. 

YY = BYX(LL) ** 2. * SXe 

XH OBKY Gas ies 25 * SY! 

SYX(LL) = (PAIRS(LL) - 1.) / (PAIRS(LL) - Deeg GSiie = ae) 
SXY(LL) = (PAIRS(LL) - i.) 7 CPAIRS(LL) = D) x) CSK2 -= 229) 
R(LL) = SUMMUL(LL) / (SUMSQ2( 1) * SUMSQ2( 2)) «** -5 
RSQ(LL) = R(LL) ** 2. 

STESTY(LL) = STDDVX(LL) * Cu, =/RSOCGLL)) *# -5 

STESTX(LL) = STDDVY(LL) + (1. - RSQ(LL)) ** 295 

Dey (it, =] WOOkK STESTY(LL) / AMEANX(LL) 

ERRDKY(LL) = DXY(LL) / (2. * PAIRS(LL)) ** .5 

DYX(LL) = 100. * STESTX(LL) / AMEANY(LL) 


BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 
BVYXY 


21 


62 
63 
64 


65 
66 
67 
68 
69 
70 
zu 
72 
73 
74 
iD 


76 
al 
78 
79 
80 
81 
82 
83 
84 
85 
86 


87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 


BVYXY100 
BVYXY101 
BVYXY102 
BVYXY103 
BVYXY104 
BVYXY105 
BVYXY106 


ERRDYX(LL) = DYX(LL) / (2. * PAIRS(LL)) ** 5 
GONBYX(LL) = SYX(LL) ** .5 / (STDDVX(LL) * ((PAIRS(LL) - Te) eee 5) EV ERE LOZ 
CONBXY(LL) = SXY(LL) ** .5 / (STDDVY(LL) * ((PAIRS(LL) - 1. )ex-5))BVYXY108 
CONAY (LL) SwveK(alsliy) ue eye a/- (PAIRS(LL) +** -5) 


at) 


CONAX (LL) ey (ial, ee ve De /, (PAIRS(LL) ** 25) 
XMAX(LL) = SPMAX(J1) * BXY(LL) + AX(LL) 


BVYXY109 
BVYXY110 
BVYXY111 


22 


48 


oS) 


50 


15 
7, 


16 
Da 


aL 
45 
46 


36 
42 
Ly 
43 
39 
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YMAX2(LL) = SPMAX(J2) * BYX(LL) + AY(LL) 


IF (LGT)49,48,49 

VDY = SPMAX(J1) / 49. 

VDX = SPMAX(J2) / 35. 

GO TO 50 

VDY = (SPMAX(J1) - SLGMIN(J1)) / 49. 

VDX = (SPMAX(J2) - SLGMIN(J2)) / 35. 
TOP2(LL) = (YMAX2(LL) - SPMAX(J1)) / VDY 
BASE2(LL) = AY(LL) / VDY 

TOP3(LL) = (XMAX(LL) - SPMAX(J2)) / VDX 
BASE3(LL) = AX(LL) / VDX 


he EGR) 5, V6, 15 
DOP? i, K 
NOt) =. TK) 
MEha eye = aC.) 
SPMAX(J1) = TSPMAX 
SPMAX(J2) = TSPMXX 
CONTINUE 

NN = NBIVAR 


IF (M+MM)36, 31, 36 

IF (81-K)45,46,46 

WRITE OUTPUT TAPE 6,19 

WRITE OUTPUT TAPE 6,18 

GO TO 39 

IF (M)43,42,43 

IF (50-K)44,46,46 

WRITE OUTPUT TAPE 6,19 

WRITE OUTPUT TAPE 6,40, (LGCOL(J) ,J=1, IM) 


BVYXY112 


BV YX YARES 
BVYXY114 
BVYXY115 
BVYXY116 
BVYXY117 
BVYXY118 
BVYXY119 
BVYXY120 
BVYXY121 
BVYXY122 


BVYXY123 
BVYXY124 
BVYXY125 
BVYXY126 
BVYXY127 
BVYXY128 
BVYXY129 
BVYXY130 


BVYXY131 
BVYXY132 
BVYXY133 
BVYXY134 
BVYXY135 
BVYXY136 
BVYXY137 
BVYXY138 
BVYXY139 


WRITE OUTPUT TAPE 6,FMT8,(R(J),J=2,NN), (BYX(J),J=2,NN), (AY(J) , J=2, BVYXY140 
1NN) , (BXY(J) ,J=2, NN), (AX(J),J=2, NN), (STESTY(J) , J=2,NN), (STESTX(J) , JBVYXY141 
2=2,NN),(DYX(J),J=2,NN) , (ERRDYX(J) ,J=2, NN), (DXY(J) , J=2,NN) , (ERRDXY(BVYXY142 
3J),J=2,NN) , (CONBYX(J),J=2,NN) , (CONAY(J) , J=2,NN) , (CONBXY(J) , J=2, NN) BVYXY143 
4, (CONAX(J),J=2,NN), (YMAX2(J),J=2, NN), (XMAX(J) , J=2,NN) ,(TOP2(J) , J=2BVYXY144 


5,NN), (BASE2(J),J=2,NN), (TOP3(J) ,J=2,NN) , (BASE3(J) , J=2, NN) 


DO 53 J=2,NN 
R(J) = 0. 
BYX(J)! ="0), 
ANICH)) = Os 
BXY(J) = 0. 
AX(J) = 0. 
STESTY(J) 
STESTX (CS) 
DYX(J) = 0. 
ERRDYX(J) = O. 
DKY(S)- =)0. 
ERRDXY(J) = 0. 
CONBYX(J) = 0. 
CONAY(J) = 0. 


O. 
OVS 


BVYXY145 


BVYXY146 
BVYXY147 
BVYXY148 
BVYXY149 
BVYXY150 
BVYXY151 
BVYXY152 
BVYXY153 
BVYXY154 
BVYXY155 
BVYXY156 
BVYXY157 
BVYXY158 
BVYXY159 
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CONBXY(J) = 0. BVYXY160 
CONAX(J) = 0. BVYXY161 
YMAX2(J) = 0. BVYXY162 
YMAX(J) = 0. BVYXY163 
TOP2( I) =1O% BVYXY164 
BASE2(J) = 0. BVYXY165 
TOP3(J) = 0. BVYXY166 
53 BASEZ(J) = 0. BVYXY167 

Cc 
IF (MM)37, 38,37 BVYSY168 
37 M=1 BVYXY169 
MM = O BVYXY170 
GO TO 32 BVYXY171 
38 NPTMP2 = NPTMP2 + 1 BVYXY172 
NROW1 = NROW2 + 1 BVYXY173 
NROW2 = NROW2 + NBIVAR - 1 BVYXY174 
GO TO 22 BVYXY175 

CG 
20 WRITE OUTPUT TAPE 6,21 BVYXY176 
CALL Hert BVYXY177 

C 
21 FORMAT(1H1,23HANOMALOUS VALUE IN DATA) BVYXY178 
18 FORMAT (1H-,107HBIVARIATE ANALYSIS (REGRESSION OF Y ON X AND X ON BVYXY179 
1Y) OF Y/X. RESULTS NOT VALID FOR COLUMNS OF LOG Y/LOG X.) BVYXY180 
19 FORMAT (1H1) BVYXY181 
4O FORMAT (1H-,93HBIVARIATE ANALYSIS (REGRESSION OF Y ON X AND X ON YBVYXY182 
1) OF LOG Y/LOG X, VALID ONLY FOR COLUMN(S),714/1H ,1314) BVYXY183 

Cc 
22 RETURN BVYXY184 

END 


Input Data Preparation for DASAN 


The following section describes the preparation of control 
and data cards and the order in which they are to be placed 
behind the program source deck. The control cards are de- 
scribed in two groups. Group I is placed directly behind the 
program source deck and is not repeated for each data deck 
included in the computer run. Group II follows Group I and 
must be repeated before each data deck in the run. This 
means that operations controlled by Group I apply to all 
data decks in a run, but that operations controlled by each 
Group II apply only to the data deck immediately follow- 
ing that appearance of Group II. 

In the description below, instructions are numbered con- 
secutively. Numbers in parentheses indicate the minimum 
and maximum numbers of control cards which are permit- 
ted for each instruction. 

All numbers punched on control cards are without deci- 
mal points, unless otherwise indicated, and are right justified 
within their fields. 


The schematic problem illustrated in Figure 1 is used as an 
example for some of the controls. All instructions are fur- 
ther demonstrated by the listing of controls in the input for 
a real problem. 


Group I Control Cards 


Fe alt) The number of data decks (from 1 to 999) 
which are to be processed in a single run. 
Punch in card columns 1-3. 

Zeal) The number of variables in the univariate 


input array (from 2 to 53), including the col- 


umn of specimen identification numbers. 
Punch in columns 2-3. 

The number is 6 for the problem in Fig- 
ure 1. 
The number of variables in the univariate out- 


put array (from 2-53), including the column 


24 
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of specimen identification numbers. (Note 
that if the number for Instruction 2 is max1- 
mal, i.e., 53, there is no room remaining for 
additional univariate variables, and the uni- 
variate output array will consist only of the 
array A shown in Figure 1.) Punch in col- 
umns 2-3. 

The number is /3 for the problem in Fig- 
ure 1. 
The number of variables in the bivariate ar- 
ray (from 0 to 56), including the column of 
specimen identification numbers. Punch in 
columns 2-3. 

The number is 9 for the problem in Fig- 
ure |. 
Signal for the calling of Subroutine RAW- 
TAB (one) or the bypassing of RAWTAB 
(zero) , punched in column 3. RAWTAB may 
be bypassed if no new univariate variables 
are to be added to the input array. 
Signal for the calling (one) or bypassing 
(zero) of Subroutine XRATIO, punched in 
column 3. XRATIO may be bypassed if no 
ratios are to be computed and if no bivariate 
analyses are to be run. 
(a) Columns 2-3: The number of univariate 
input variables (from 0 to 30) from which a 
constant, the value of which is specified in In- 
struction 33, is to be subtracted. 
(b) Columns 5-6: The number of univari- 
ate input variables (from 0 to 30) which are 
always to be divided by a constant, the value 
of which is specified in Instruction 34. (See 
also Instruction 9.) 
(c) Column 9: The number of univariate 
input variables (from 0 to 5) which may or 
may not be divided by the constant referred 
to above (7b) depending on the value of a 
signal included with the data for each speci- 
men. (See discussion of data deck in Instruc- 
tions 21 and 37.) 
(d) Columns 11-12: The number of uni- 
variate variables (from 0 to 30) to be gener- 
ated by subtracting one input variable from 
another. 
(e) Columns 14-15; The number of uni- 
variate variables (from 0 to 30) to be gener- 
ated by adding two input variables. _ 
(f) Columns 17-18: The number of uni- 
variate variables (from 0 to 30) to be trans- 
formed to logarithms to the base 10 and to be 
added as new columns to the original array. 


8. 


10. 


1 


(0-2) : 


(0-2): 


(0-1): 


(0-2): 


(0-2) : 


(0-2) : 


(0-2) : 


For the problem in Figure 1, the control card 
for Instruction 7 is punched 01 01 00201 02. 
The column numbers of variables in the 
univariate output array from which a con- 
stant, specified in Instruction 33, is to be sub- 
tracted. (See also Instruction 7a.) Punch 
numbers in fields of 3 in columns 1—72. 

The card is punched 006 for the problem 

in Figure 1. 
The column numbers of variables in the 
univariate output array which are always to 
be divided by a constant, specified in Instruc- 
tion 34. (See also Instruction 7b.) Punch 
numbers in fields of 3 in columns 1-72. 

The card is punched 004 for the problem 

in Figure 1. 
The column numbers of variables in the um- 
variate output array which may or may not 
be divided by a constant, specified in Instruc- 
tion 34, depending on a signal read in with 
the data for each specimen (see Instruction 
21). The maximum number of column num- 
bers which can be listed here is 5. Punch 
numbers in fields of 3 in columns 1-15. 

No card is required here for the problem in 

Figure 1. 
Column numbers of variables in the um- 
variate output array which are to serve as 
minuends (see Instructions 7d and 12). Punch 
in fields of 3 in columns 1-72. 

The card is punched 002002 for the problem 

in Figure 1. 
Column numbers of variables in the univari- 
ate output array which are to serve as sub- 
trahends, so that the first subtrahend specified 
here is subtracted from the first minuend 
specified in Instruction 11, the second subtra- 
hend from the second minuend, etc. Punch 
numbers in fields of 3 in columns 1-72. 

The card is punched 002002 for the problem 

in Figure 1. 
Column numbers of variables in the univari- 
ate output array to which other variables, 
specified in Instruction 14, are to be added. 
(See also Instruction 7e.) Punch in fields of 3 
in columns 1-72. 

The card is punched 005 for the problem in 
Figure 1. 

Column numbers of variables in the univari- 
ate output array which are to be added to 
variables specified in Instruction 13, such that 
the first variable specified here is added to the 


15. 


16. 


liye 


(0-2): 
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first variable specified in Instruction 13, the 
second variable to the second variable, etc. 
Punch in fields of 3 in columns 1-72. 

The card is punched 008 for the problem in 
Figure 1. 
Column numbers of variables in the univar- 
late output array which are to be transformed 
to logarithms to the base 10 and which are to 
be entered as new columns in the univariate 
output array. (See Instruction 7f.) Punch in 
fields of 3 in columns 1-72. 

The card is punched 002003 for the prob- 
lem in Figure 1. 
(a) Columns 2-3: The number of univariate 
variables which are to be fitted on an output 
page. (This includes a column of specimen 
identification numbers on each page.) 
(b) Columns 5-6: The number of segments 
into which the univariate output array must 
be broken in order to fit on standard output 
paper. Each output segment of the array be- 
gins with a column occupied by the specimen 
identification numbers, so that rows in the 
segment may be readily identified. Calculate 
the number of segments by means of the fol- 
lowing formula: 

LG, 

where N, the number of univariate output 
segments, is rounded to the next highest whole 
number; 7 is the number of variables in the 
univariate output array (Instruction 3) ; and 
C is the number of columns which can be 
fitted on an output page (Instruction 16a). 

For the sample output (see listing), 7" is 
6, and C is 14, so that N is 3, after rounding 
2.7 to the next highest whole number. 
(a) Columns 2-3: The number of bivariate 
variables which are to be fitted on an output 
page. (This includes a column of specimen 
identification numbers. ) 
(b) Columns 5-6: The number of segments 
into which the bivariate array must be broken 
in order to fit on standard output paper. Cal- 
culate by means of the formula in Instruction 
16b, where T is the total number of variables 
in the bivariate array (Instruction 4) and C 
is the number of bivariate columns which are 
to be fitted on an output page (Instruction 


Mia) 


18. (0-10) : 


19: 


(2n): 


For the sample output (see listing) , J’ is 42 
and C is 12, so that N is 4 after rounding 
3.8 to the next highest whole number. 
Indexing of column numbers of variables to 
be used in computation of ratios, y/x, begin- 
ning with the first ratio to appear on the 
first bivariate output page and proceeding 
from left to right through consecutive pages. 
Each ratio is represented by a group of 12 
card columns, punched as follows: 

(a) Columns 2-3: Column number of y in the 
univariate output array. 

(b) Columns 5-6: Column number of x in the 
univariate output array. 

(c) Columns 8-9: Column number which the 
ratio will have on its respective output page. 
(The column of specimen identification num- 
bers is column | on each page. ) 

(d) Column 12: Signal for the log (base 10) 
transformation of both y and x in order to 
produce bivariate statistics of the ratio both 
before and after log transformation. Punch 
zero if no transformation; punch one for 
transformation. 

Indexing for the next ratio begins in the 
next field of 3 columns, and so on, through 
the first 72 columns of the card (6 ratios in- 
dexed per card). Indexing then begins again 
on a second card. Because the maximum num- 
ber of ratios permissible is 55, the maximum 
number of cards allowed here is 10. 

If no ratios are to be computed (i.e., if the 
control number described in Instruction 4 is 
zero), no cards should be present for this in- 
struction. 

For the schematic problem in Figure 1, two 
cards would be required which would be 
punched as follows: first card: 02 03 02 01 
03 04 03 01 06 03 04 00 06 05 05 01 09 04 06 
01 10 04 07 00; second card: 02 06 08 01 12 
13 09 00. 

Names of variables to be placed as column 
headings on each output page. The column 
headings of each page occupy two cards (col- 
umns 1-72 of the first card and columns 1-48 
of the second card), with the spacing begin- 
ning with card column 1, corresponding to 
print positions on the output page (120 print 
positions per line). The names of variables on 
consecutive pages occupy consecutive pairs of 
cards. The number of cards required here is 


26 


22. 


(1-5): 
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twice the number of output pages. These cards 
must be present even if they are blank. 

Numbers of cards to be read in for each of 
8 variable formats 2128) 
Punch in fields of 3 in card columns 1—24. 
The number of cards for each format must be 


(Instructions 


at least one and cannot exceed 5, with the 
exception of the eighth format, for which up 
to 8 cards may be present. 

The format, in FORTRAN II, for the read- 
ing in of data cards. The format must begin 
with a left-parenthesis and end with a right- 
parenthesis. 

If the control described in Instruction 7c 
is greater than zero, the following applies: 

Make provision at the beginning of each 
data card (or each group of cards if a speci- 
men is represented by more than one card) 
for reading in the signals which will determine 
whether certain variables (the column num- 
bers of which were specified in Instruction 10) 
are to be divided (—1) or not divided (+1) 
by a constant, the value of which is specified 
in Instruction 34. 

There can be no more than 5 signals, and 
the signals are to be read in as fixed-point 
(integer) variables by means of an “I” specifi- 
cation as shown in the format listed with the 
controls for the sample output. All other 
variables on the data card must be read in 
as floating-point variables by means of an 
ae LegnOrmaeb ye 

In addition to the example given by the 
listing of control cards for the sample input, 
the following is given: 


specification. 


A data card which is punched, beginning in 
the first card column, + 1—1—1001760321947, 
and which is read in by the format (312, 
F3.0, 2F3.1, F3.3), will be read in as the 
following values: +1, -1, Pl. 76.0) 32510 
and .947. 

For further instruction on format specifi- 
cation in FORTRAN II, particularly with 
regard to Hollerith “H” specifications and 
carriage control, the reader is referred to the 
FORTRAN General Information Manual 
(IBM, 1961) or to McCracken (1961 Nie 
The format, in FORTRAN II, for the 
punching out of the univariate output array. 
The order in which the variables are punched 


oo 


On 


26. 


Dies 


28. 


(1-5): 


(1-5) : 


. (1-5): 


(1-5): 


(1-5): 


out is the same as the order of the subarrays 
shown in Figure 1. 

This is also the format for the reading in of 
maxima for the variables in the univariate 
output array (Instruction 29) and for the 
reading in of log minima for the variables 
which will be transformed to logs by Instruc- 
tion 18d. 

The format, in FORTRAN IT, for punching 
out the bivariate array. 

The format, in FORTRAN II, for printing 
the univariate output array. The format is 
repeatedly used for each row of the array 
throughout all segments, so that the maximum 
number of decimal places which will be 
required by any one variable should be used 
for all variables. The same applies to Instruc- 
tion 25. 

The format, in FORTRAN II, for the print- 
ing of the bivariate array. (See Instruction 
245) 

The format, in FORTRAN II, for printing 
the univariate statistics beneath each table of 
the univariate array. Column alignment is 
achieved by matching the field specifications 
to those of Instruction 24. Note in the listing 
of the sample input that an “I” (integer) 
specification must be used for the row of 
numbers of non-zero values. 

The format, in FORTRAN II, for printing 
the univariate statistics beneath each table of 
the bivariate array, again with “IT” specifica- 
tions for the row of numbers of non-zero 
values. 

The format, in FORTRAN II, for printing 
the bivariate statistics beneath each segment 
of the bivariate array. 

Card(s) containing the anticipated maximum 
values of all the variables in the univariate 
output array of this run or any future run 
analyzing the same variables. The maxima 
are hypothetical and are specified in order to 
determine the calibration of the axes of the 
machine-plotted bivariate scatter diagrams, 
which in turn allows the computation of 
slopes and intercepts in terms of the scale- 
units of the diagrams. (See Instructions 11 
and 12 in the description of input data for 
VPLOT.) In this manner, the calibration is 
kept uniform throughout all runs so that the 


30. 


Bile 


(lea): 


(1): 
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scatter diagrams for any given ratio can be 
superimposed for comparison. 

The axes of plots of non-log data are as- 
sumed to cross at (0,0). For log data, how- 
ever, the axes may cross at (a,b), where a and 
b are minima specified for x and y on control 
cards (Instruction 30). 

The punching format, and hence the num- 

ber of cards occupied by the maxima, is de- 
termined by Instruction 22. There must be at 
least one card present, even if no maxima 
are specified. 
Card(s) containing the anticipated minimum 
log values of those variables which will be 
transformed to logarithms (base 10) by In- 
struction 18d. These minima are specified in 
order to determine the coordinates of the end 
points of regression lines in terms of scale 
units of a scatter diagram plotted by VPLOT. 
Generally the specified log minima are -l. 
or —2. 

The punching format is the same as in 
Instruction 29, with no data required in the 
fields of those variables which will not be 
transformed to logs. 


Group II Control Cards 


(a) Columns 1-2: Signal for the punching 
out of the univariate and bivariate arrays, as 
follows: 
00=no punch-out 
+1=punch univariate array only 
+2=punch bivariate array only 
+3=punch both univariate and bivari- 
ate arrays. 
(b) Columns 4-5: Signal for the printing out 
of the univariate array and accompanying uni- 
variate statistics with the order of the speci- 
mens in the data deck preserved: 
00=no print-out 
+1=print. 
(c) Columns 7-8: Signal for the printing 
out of the bivariate array with its accompany- 
ing univariate and bivariate statistics, with the 
order of specimens in the data deck preserved : 
00=no print-out 
+1=print 
(d) Columns 10-11: Signal for the printing 
out of the univariate array and accompanying 
univariate statistics, with the specimens sorted 


2 
oO 





32. 


335 


oa: 


3D; 


36. 


Si 


(1): 


(1): 


()e 


according to decreasing values of a univariate 
variable specified in Instruction 31f: 

00=no print-out 

-+1=print. 
(e) Columns 13-14; Signal for the printing 
out of the bivariate array and accompanying 
univariate and bivariate statistics, with the 
specimens sorted according to decreasing val- 
ues of a univariate variable specified in In- 
struction 31f. 

00=no print-out 

+1=print 
(f) Columns 16-17: Column number of a 
variable in the univariate output array, the 
decreasing values of which will control the 
ordering of specimens (rows) in the output, 
if such output is indicated by the options de- 
scribed in Instructions 31d and 3le. 
Any name or sample-identification informa- 
tion to be printed at the top of each output 
page. Punch in columns 1-72 of the first card 
and columns 1-48 of the second card. 
A constant which is to be subtracted from 
certain univariate variables specified in In- 
struction 8. Punch in columns 1-6 without 
decimal point. The decimal will be placed by 
the computer before the last two digits; for 
example, 017291 will be read in as 172.91. 
A constant by which certain univariate varl- 
ables, specified in Instructions 7b and 9, are 
to be divided. Punch in columns 1-6 in the 
same manner as the constant described in 
Instruction 33. 
A signal allowing for the following choices of 
lines to be printed at the top of the first 
output page, below the sample-identification 
information : 

—1: x LEFT VALVES 

00: x SPECIMENS 

+1: x RIGHT VALVES 
where x is the number of specimens in the 
data deck. Punch signal in columns 1-2. 
The number of specimens (from 1-100) in 
the data deck. Punch in columns 1-3. 


Data Deck 


The data cards must be punched according 
to the format specified in Instruction 20. 
There can be any number of continuation 
cards for each specimen provided that these 
are specified by the format. 
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Alternative Statements and Subroutines for DASAN 


Comments following statement No. 28 in Subroutine 
RAWTAB mark the position in which may be inserted 
FORTRAN II statements for the computation of special 
variables. As shown in Figure 1, these variables form a sub- 
array B within the univariate output array. One special- 
computation statement and the two statements that must 
accompany it are included in RAWTAB for the purpose of 
illustration. (The “C” 


in column 1 makes each of these 


statements inoperative in the program.) In this case a vari- 
able named “P.AR” in the sample output is computed, where 
P.AR is the area of a triangle, the base and height of which 
appear in the univariate output array in columns 21 and 20 
respectively. No special control cards are required, but the 
number of variables in the univariate output array must be 
increased by / for each special computation carried out and 
adjustments must be made in the output formats. 


Sample Input Data for DASAN 


Listed below is the entire assembly of data for a test run of 
DASAN. The problem involves the analysis of 12 specimens, 
with 28 variables in the univariate input array, 36 variables 
in the univariate output array, and 42 variables in the 
bivariate array. The numbers on the far right side of each 
card refer to instruction numbers in the discussion of input 


preparation. Included within the listing (Instruction 28) 
is the format for the print-out of bivariate statistics involv- 
ing regressions. The alternative format, to be used with the 
bivariate subroutine involving the reduced major axis, 1s 
shown below the input list. 


DATA 
001 a 
028 2 
036 B 
o42 4 
001 5 
71 6 
0100 5 02 02 03 7 
012 fs) 
013018019020021 10 
006007 12 
005011 We 
004008 13 
002007 14 
029012006 1 
014003 16 
012004 17 
06. 32: 02 0006/28 103 00 04 O2 OF 00 05 30 05 00 07 02 06 00 07) 52 07,00 18 
08 04 08 00 08 32 09 00 07 08 10 00 33 32 11 00 09 06 12 OC 10 06 O02 00 
31 O07 03 00,14 32.04 00 12 06 05 01 12 32 06 00) 13 32 07 00 13 0608 00 
18) 19) 09.00 18) 06 10 00 20 06 11 00 21 52 12 00 20 21°02 00 22 32,105,100 
23.06 O04 00: 22 25 05 00,24 06 06 00 25 32-07 00 24 25 08 00 26 06 09°00 
27 32 10 00 05 03 11 00 05 O4 12 00 31 02 02 00 34 36 03 00 29 06 O4 00 
07 06 05 00 12 28 06 00 O02 03 07 00 29 12 08 00 07 28 09 00 
SPEC DG GP AD AK AM DF CD 19 
EI BJ DE LO RIB HT HI 
SPEC PLICAE ANTCOS POSCOS INTRIB INTGRV (AD) (CE) ( 
IK ) (PQ) (BJ) (AB) (MN) (FG) 
SPEC ALM P. AREA KM EF AG CF LOG P.AR LO 
GLO LOG AM 
SPEC AM/AG AM/ALM AD/DG AK /KM DF/DG DF/AG CD/A 
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D CD/AG DF/CD CF/AG EI/AM 
SPEC BJ/AM EF/DF HI/AG LO/AM LO/AG RIBHT/AG RIBHT/ 
AM  RIB/GRV- RIB/AM (AD) /AM (CE) /AG 
SPEC (AD)/(CE) (IK)/AG (PQ)/AM (IK)/(PQ) (BJ)/AM (AB)/AG (BJ)/( 


AB) (MN)/AM  (FG)/AG AK/GP AK/AD 
SPEC EF/DG LOG PA/AM PA/AM DF/AM LO/ALM DG/GP P.AR/L 
0 DF/ALM 


001001001001001003003008 
(Sie hs Ones .d ol OFS... 5F210, 485.1 /6F 32, Fe. 1) 
(HE O.FG.1, 2Or5) 1, 9X, VHA/I4E5.1,9X, LHB/5F5. 1, 5F7 +, 33X, LHC) 
(F4.0,11F6.3,9X,1HA/a2F6.3,7X,1HB/12F6.3,7X, 1HC/6F6.3,43X,1HD) 
GHO/ GH Fa. O;hS 2. 12F9.1)),) 
(1HO/(1H ,F4.0,F9.3,10F10.3)) 
(1H-,4HNUM.,16,1219/1H ,4HMAX.,F8.1,12F9.1/1H ,4HMIN.,F8.1,12F9.1/1H , 5 
HRANGE,F7.1,12F9.1/1H ,4HMEAN,F8.1,12F9.1/1H 5HSTDDV,F7.1,12F9.1/1H ,5HS 
DVMN,F7.1,12F9.1/1H ,3HV ,13F9.1) 
(1H-,4HNUM.,18,10110/1H ,4HMAX.,F9.3,10F10.3/1H ,4HMIN.,F9.3,10F10.3/1H 
, DHRANGE, F8.3,10F10.3/1H ,4HMEAN,F9.3,10F10.3/1EH ,5HSTDDV,F8.3,10F10.3/ 
1H ,5HSDVMN,F8.3,10F10.3/1H ,3HV ,11F10.3) 
GH, SHR 4 222P10/53/1 HO. 5HA (YX) ,F8.3, LOF10.3/1H.., 4HBCY) «FO. 3..110F10.371.H0, 
SHAGYO hie -SeOPlOa5/1n ehHB(X)n 69 55, 1I0F1L0s 3/1 H0,6HS CY x) a7 Ss LOR1O GS /1 
Her GUSOGy) ib 7. 5 1OR1O!. 3/UHO-6HDCY .X)),F7 55. 1CRLOs3/iH », OHS GDYX)), £7. 3.10 
F1053/1H ,6HD(X.Y),F7.3,10F10.3/1H ,6HS(DXY) ,F7.3,10F10.3/1HO,6HCONAYX,F 
7.3,10F10.3/1H ,5HCONBY,F8.3,10F10.3/1H ,6HCONAXY,F7.3,10F10.3/1H ,5HCON 
BX,F8.3,10F10.3/1HO,6HYMAX 2,F7.3,10F10.3/1H ,4HXMAX,F9.3,10F10.3/1HO, 6H 
TOP VY E75. 1OFTO.5/1H ,OHBASEYX, F7.3,10F10.3/1H ,GHIOP XY,F7.3, 10F10.3 
VAL ROHBASEXY., F725, 10F1O.3) 
GOmeerOnO Ft OLO0 4550.40.00 75.0 25.0 27.0 14.50 19.5 25.20 26.0) 3.5 50-0 
ZOROMIO Ol IeHOMNS JF R550) 565 465542750 2iksOv sIe 0 AS (042020 9F 0110050 
2OROMS> 40M 65505409 96.5 1501071 41150 1.8752 
98. =1.0 -1.0 


BIL OO MOON! tl. NOG 
AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLA., MCZ 197942 ,232557,232 
553 -- LIVING 


000945 

002675 

+1 

012 

=a les S01 S150 280 80 43 61 66186 1615519 3 4560370410460 
340 

et=l=Pote e2algqimiolesi2s 240 72.71 35,61 “50N74. 1305419) 3 5390300540290 
290 

atat=V—o1) SlTOLO9IA ZI21 219-70 64.30 51) -43160) 11915214 3: 6550240280320 

19 260 

+1-1-1-1-1 4139141142151 274 89 89 42 69 59184 1515518 3 6410270400410 
355 

+1-1-1-1-1 5116116121120 239 80 78 38 63 50166 1313421 4 6300250330380 
285 


+1-1-1-1-1 6127127135141 260 82 79 42 67 56180 1414619 3 6430280380440 
515 


29 


20 


22 
oo 
24 
25) 
26 


ei 


28 


29 


30 


30 
Sil 
32 


33 
34 
35 
36 
37 


30 


tejel—et—1) 7251735129 2435 75 
24 290 

panel is Oded dT MO SIAL AIS 
270 

fil Oe 99e 96100105 
46 21 220 

edi) tl LOM Oud) 5125 
66 54 80 49 49 17 270 

#1-1-1-1-1 11127121126125 
85 66 95 49 66 320 

+P=—1—-2=) 12726111122130 


66 27 320 


(eme2 


222 70 65 34 
66 
68 
86 


81 


193 64 31 


226 66 34 
229 C2557, 


236 On o> 


(1H-,3HR ,11F10.3/1HO,3HA 


DD 
29 
50 
50 
68 
62 


slIFLO.3/1H ,6USIECA) F725, 1OFI0. 3/150, 5HB 
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50171 1413219 3 6420300350410 


51162 1312521 3 6350210260320 


41151 1211120 2 6290180270280 


49161 1212120 3 6350200290310 


57185 1414420 3 6450190370520 


60188 1614120 3 6400210410460 


28 


,11F10.3/1HO,6HD( RMA) ,F7.3,10F10.3/1HO,6HMEAN Y,F7.3,10F10.3/1H ,6HMEAN 
KX, B75, LOFTO.5/1HO0.6HYMAX. 1, F735. 1LOF10. 3/7180, 54 TOP.1. 6.5, JOFIO. 3/18) 6H 


BASE, F723. TOr1l0¢3)) 


Sample Output from DASAN 


Listed below are the punch-outs of the univariate output 
array and bivariate array and the printed output for the 
sample problem. In the printed output, variable names con- 
sist of at least two letters, with parentheses enclosing letters 
which represent lower-case letters. The abbreviations for 
statistics have already been explained in Table 1 and the 
use of the output in the plotting of reduced major axes and 
regression lines will be explained at the end of the descrip- 
tion of VPLOT. 

The specimens, which were ordered from 1 through 12 


PUNCH-OUT OF UNIVARIATE OUTPUT ARRAY 


in the input array, have been rearranged in order of de- 
creasing values of the variable AM. As in the simplified 
sample problem illustrated in Figure 1, certain columns of 
statistics are not valid. In arrays of bivariate statistics these 
columns are indicated in output messages. In the arrays of 
univariate statistics for bivariate data, the statistics for va- 
riable No. 36 (Log P.A./Log AM) are or may be incorrect. 

The bivariate analyses from an additional run using the 
alternate bivariate subroutine are shown for the second to 
the last page of the output. 


Ve ee Ss Or MO a>. e2oc On =O. SiO ees) SOc 6 s6u  Oaln coe on S 
1OFON SoOr ek Or 201 De eb ey =O On OR On =O. POR, S40 
Wag sO, NO BeO.> On OallL99 OL 96a i 47? 

2 SPA SO" VADs 2G eh On “72 Ra Seb Geil 50) Tua) © dash alsaieen 
UG ROe 5 On SOw 1s 5" 1S eS 5) O04. =O 2-0 =On One =O, FeOnO 
0.9 11.2 2.2 24.4 14.3-0.0331 0.9004 1.3802 

Sie MO LOO el e778 eee On. PLO) bic i Ol bind ole Ss iGmGn al.me alan 
ZO sO eG On eles. LORGe Ol IZ Oe sOe me O Len On aOr 19) 26.<0 
0.6 9.8 2.7 22.7 13.4-0.2034 0.8162 1.3404 

He) 9) ee Sl 27 SO OO 42629) = 5) Ou OnOue daca ali 5 
HO sOnm Se OlenG «Onley Oneulin Soule, 5a —O).8 ¥—O= Oe Olen Om a0 mao eo 
Voled2.5° 320.26. 2726. O-0592 O.95ue 12k578 
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0.226 0.324 0.516 0.315 0.299 0.055 0.058 1.400 0.065 0.054 0.060 0.854 
=O. -0. Ol -0. -0. O% -0. 0.094 1.103 05970-04195) 0.001 
O2021N02 505 O.2649 b-05is O51Sr O2255 

Bl 0.945 102822 1.117 2.075 0.651 06298 O5524 0.277 1.077 10.570 koe 
0.266. 0.271 0.532 0.30% 0.287 0.055 0.059 1.667 0.059 0.044 0051 O.812 
=O), -0. On -0. -0. Oo: -0. -0. 1.095 0.927 0.171-0.175 
0.026 0.315 0.250 1.057 0.086 0.259 

G9. 0.970 0-877) 1.010 12193 0.667 06552 0.64010. 522 1.051, 0.655, 10.1615 
0.259 0.379 0.558 0.293 0.284 0.060 0.062 1.611 0.056 0.052 0.053 0.964 
-0. -0O. Of -0. -0. 05 0.238 0.106 15071 1.4050 02255-0216 
02027 0.342 02257 1.010 070935 0.500 

10. 1.00% 0.837 1.045 1.19% 0.616 0.302 0.574 05293 1.050. 02596) O 1150 
0.221 0.279 0.538 0-294 0.296 0.053 0.053 1.750 0.058 0.046 0.052 0.955 
0.293 0.2591 222 0.3554 0.218 15653) 05217 0.076 15108 15070+0.075-0.149 
0.028 0.301 0.246 0.991 0.094 0.252 

12, 0.945 0.747 0.992 1.096 0.677 0.340 0.651 0.324 1.049 0.664 0.155 
0.285 0.337 0.569 0.379 0.358 0.055 0.059 2.368 0.070 0.058 0.077 0.712 
06.356 0.276) 1..288"0.3597.0.194 1.959)0.276-0.. 1.033 0.992 0.228 0.093 
0.056 0.360 0.283 1.050 0.149 0.269 

1), 0.952 00737 0.968 1.226 0.643. 0.527 ©:625 ©2306. 1.066 0.6535 O.048 
0.263 0.259 0.569 0.396 0.377 0.065 0.068 1.905 0.063 0.065 0.069 0.891 
-0. -0. Ox -0. -0. Or 0.280 0.1109" 1LL7D W066 102167) 0.087 
O2056 O. 343 Os292) 1.155 -O.T4) 05253 


AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLAs, MCZ 197942,232557,232553 -- LIVING 


12 RIGHT VALVES 


SPEC 0G GP AD AK AM DF co EI BJ OE Lo RIB HT 
Le 14.5 13.0 14.8 15.0 28.0 -0. 8.0 4.3 661 6.6 9el 1.6 
Ge 13.9 14.1 14.2 15.1 27.4 8.9 8.9 4.2 6.9 529 8.9 1.5 
6. 12.7 12.7 13.5 14.1 26.0 8.2 7.9 4-2 6.7 526 8.6 1.4 
Te 12.3 11.7 13.3 12.9 24.3 1.4 7.2 3.5 5.5 5.0 7.6 1.4 
2. 11.9 11.0 12.5 12.8 24.0 7.2 Tel 3.5 661 5.0 7.9 1.3 

ll. 12.7 12.1 12.6 12.5 23.9 8.6 8.2 3.7 6.8 5.7 9-1 1.4 
5. 11.6 11.6 12.1 12.0 23.9 8.0 7.8 3.8 6.3 5.0 Tel 1.3 
L2. 12.6 ll. 12.2 13.0 2326 8.1 726 3.5 6.2 6.0 924 L.6 
10. 11.0 Aled 11.5 12.3 22.6 6.8 6.6 304% 5.0 409 6.6 1.2 
8. 11.1 10.5 12.4 11.5 22.2 Ted 6.5 34 5.9 Sel 6.7 1.3 
36 11.0 10.9 11.7 12.1 21.9 7.0 6.4 3.0 5el 43 6.6 tel 
9 9.9 9.8 10.0 10.5 19.3 6.6 6.4% 3-1 5.0 4el 5.6 1.2 

NUM. 12 12 12 12 12 LL 12 12 12 12 12 12 

MAX. 14.5 14.1 14.8 15.1 28.0 8.9 8.9 4.3 6.9 6.6 9.4 1.6 

MIN. 9.9 9.8 10.0 10.5 19.3 6.6 6.4% 3.0 5.0 4el 5.6 1el 

RANGE 4-6 4.3 4.8 4.6 8.7 23 2.5 1.3 1.9 2.5 3.7 0.> 

MEAN 12.1 11.6 12.6 12.8 23.9 7.6 7.4 326 6.0 5.3 7.8 1.4 

STOOV 1.3 1.2 1.3 1.4% 224 0.8 0.8 0.4 0.7 0.7 1.2 Ved 

SOVMN 0.4 0.3 0.4% 064 0.7 0.2 0.2 O.l 0.2 0.2 0.4 UU 

v 10.7 10.2 10.2 10.6 10.0 10.2 11.0 11.6 11.5 13.7 15.9 11.5 
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AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, 


12 RIGHT VALVES 


SPEC 


1. 
4. 
6. 
Te 
2. 
ll. 
5. 
12. 
10. 
8. 
3. 
9. 


NUM. 
MAX. 
MIN. 


RANGE 


MEAN 


SITDOV 
SDVMN 
v 


AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, 


12 RIGHT VALVES 


SPEC 


1. 
4. 
66 
Te 
26 
ll. 
5. 
l2. 
LO. 
8. 
36 
Ge 


NUM. 
MAX. 
MIN. 


RANGE 
McAN 

STDOV 
SOVMN 


Vv 


PLICAE 


19.0 
18.0 
19.0 
19.0 
19.0 
20.0 
21.0 
20.0 
20.0 
21.0 
21.0 
20.0 


12 
21.0 
18.0 
320 
19.7 
1.0 
0.3 
409 


ALM 


34.0 
33.5 
31.5 
29.0 
29.0 
32.0 
28.5 
32.0 
27.0 
27.0 
26.0 
22.0 


12 
34.0 
22.0 
12.0 
29.3 
3.5 
1.0 
11.9 
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11.9 
10.6 
10.3 
10.7 

9.8 

8.8 


13.0 
8.8 
42 

11.1 
1.2 
0.3 

10.4 


RRR RR RR RR REN 
Chielislehetelnoiemehelience 


~- 
TOOK KK NN 
mH WewWURE NUON eR 


oeeeeee 
Ce wvorr 


~ 


m 
nn 


NNR RF NUWUNNNNWO 
eee eee ee © oO Oo 


- 
SOONK Kee 

eoeeeeeve 

CeFrur OO WN DOORKOCONFSAO 


~ 


FLAs, MCZ 197942, 232557,232553 -- LIVING 
INTGRY (Ad) (PQ) 
1.4 1.5 1.7 -0. 
1.0 1.5 1.5 -0. 
1.0 1.4 1.6 -0O. 
lel 1.3 1.5 -0O. 
1el 1.3 1.5 -0O. 
0.7 1.4 1.9 6.6 
0.9 1.2 1.4 -0. 
0.8 1.5 1.7 -O. 
0.7 lel 1.2 564 
0.8 1.0 1.2 -0. 
0.9 1.) 1.2 =Oie 
0.7 1.0 1.0 -O. 
12 12 12 2 
1.4 1.5 1.9 6.6 
0.7 1.9 1.0 5.4 
0.7 0.6 0.9 1.2 
0.9 1.3 1.5 6.0 
0.2 0.2 0.3 0.8 
OL O.l O.1 0.6 
22.6 16.2 18.6 1 14.1 
FLAs, MCZ 197942,232557,232553 -- LIVING 
AG CF LOG AM 
29.3 O. 0.1 1.0 1.4 
28.1 17.8 O.1 1.0 1.4 
26.2 loel 0.1 0.9 1.4 
25.6 14.6 0.0 0.9 1.4 
2464 14.3 -0.0 0.9 1.4 
25.3 16.8 O.L 1.0 1.4 
23.7 15.8 -0.1 0.9 1.4 
24.8 15.7 0.1 1.0 1.4 
22.5 13.4 -0.2 0.8 1.4 
23.5 13.5 -0.2 0.8 1.3 
22.7 13.4 -0.2 0.8 1.3 
19.9 13.0 -0.3 0.8 1.3 
12 lL 12 12 
29.3 17.8 1.0 1.4 
19.9 13.0 0.8 1.3 
924 4.8 0.2 0.2 
24.7 14.9 0.9 1.4 
2.5 1.6 0.1 0.0 
0.7 0.5 0.0 0.0 
10.3 10.7 8.0 3-2 
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AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLA«, MCZ 197942,2325574232553 -- LIVING 


12 RIGHT VALVES 


SPEC AM/AG AM/ALM AOD/0G AK/KM OF /0G DF /AG CD/AD CO/AG OF/CD CF/AG EI/A™M 
le 0.956 0.824 1.021 1.154 -O. -O.~ 0.541 0.273 -0. Oo. 9.154 
4 0.975 0.818 1.022 1.228 0.640 0.317 0.627 0.317 1.000 0.633 Jel53 
6. 0.992 0.825 1.063 1.185 0.646 0.313 0.585 0.302 1.038 0.615 9-162 
Tl. 0.949 0.838 1.081 1.132 0.602 0.289 0.541 0.281 1.028 0.570 Jel 44 
2. 0.98% 0.828 1.050 1.143 0.605 0.295 0.568 0.291 1.014 0.586 Uel46 

tl. 0.945 0.747 0.992 1.096 0.677 0.340 be651 0.324 1.049 0.664% 3-155 
Se 1.008 0.839 1.043 1.008 0.690 0.338 0.645 0.329 1.026 0.667 Uel59 
12. 0.952 0.737 0.968 1.226 0.643 0.327 0.623 0.306 1.066 0.633 0.148 
LO. 1.004 0.837 1.045 1.194 0.618 0.302 0.574 0.293 1.030 0.596 92152 
8. 0.945 0.822 Le1lL7 1.075 0.631 0.298 06524 O.277 1.077 0.574 9.153 
3e 0.965 0.842 1.364 1.235 0.636 0.308 0.547 0.282 1.094 0.590 Oel37 
Je 0.970 0.877 1.010 1.193 0.667 0.332 0.640 0.322 1.031 0.653 C.161 

NUM. 12 12 12 12 il ll 12 12 Ll ll 12 

MAX. 1.008 0.877 Lell7 1.235 0.690 0.340 0.651 0.329 1.094 0.667 Uel62 

MIN. 9-945 0.737 0.968 1.008 0.602 0.289 0.524 0.273 1.000 0.570 OeL3T 

RANGE 0.06% 0-140 0.149 0.226 0.0868 0.051 0.127 0.056 0.094% 0.096 U2025 

MEAN 0.970 0.820 1.040 1.156 0.641 0.314 0.589 0.300 1.041 0.617 Je152 

STODV 0.023 0.039 0.040 0.069 0.028 0.018 2046 0.020 0.028 0.036 0-0U/ 

SOVMN 0.007 0.011 0.012 0.020 0.008 0.005 0.913 0.006 0.008 0.011 0.007 

v 2.338 4.796 3.883 52986 4.370 5-614 7.826 6.594 2-680 5.768 42143 


BIVARIATE ANALYSIS (REDUCED MAJOR AXIS) OF Y/X- RESULTS NOT VALID FOR COLUMNS OF LOG Y/LOGX. 


K 0.974 0.916 0.934 0.832 0.905 0.836 0.730 0.806 0.975 0.844% 0.922 
A 0.948 0.686 0.983 1.182 0.700 0.357 0.637 0.321 0.936 0.734 Je176 
SIG(A) 02062 0.079 0.102 0.189 0.090 0.059 0.126 0.055 0.063 0.119 9.02° 
B 0.530 3-821 0.675 -0.310 -0.700 -1.045 -02623 -0.539 0.757 -2.857 -0.57d 
D(RMA) 2.297 4.597 3.804 6.091 4.167 5.191 7.637 6.441 2.401 5.276 32982 
MEAN Y 23.925 23.925 12.567 12.817 7.618 7.618 7.383 7.383 7.618 14.945 3.633 
McAN XK 24.667 29.292 12.100 11.108 11.882 242245 12.567 24.667 71.327 242245 23.925 
YMAX 1 81.2149 726454 39.987 41.049 27.302 29.326 28.047 26.760 26.039 59.554 12.624 
TOP 1 4.017 -1.664 ie 49) 1.286 4.513 8.478 1.900 -0.435 2.037 414.201 4.814 


BASE Ll J.346 2.497 0.735 -0.380 -1.371 -2.048 -1.131 -0.978 1.484 -22.219 -2.024 


AEQUIPECTEN NUCLEUS -- MIAMI 


12 RIGHT VALVES 


SPEC 


le 
a. 
6. 
Te 
2. 
Al. 
Se 
12. 
10. 
8. 
3. 
9. 


NUM. 
MAX. 
MIN. 
RANGE 
MEAN 
STODV 
SOVMN 
v 


BIVARIATE ANALYSIS 


A 
SIG(A) 


8 


D(RMA) 


MEAN Y 
MEAN X 


YMAX 1 


rop i 
BASE 1 


BJ/AM 


0.218 
0.252 
0.258 
0.226 
0.254 
0.285 
0.26% 
9.263 
0.221 
3-266 
0.233 
0.259 


12 
0.285 
0.218 
0.067 
0.250 
0.021 
0.006 
8.241 


0.701 


0.285 
92059 


-0.847 
7.839 


5-967 
23.925 


20-512 


2-542 
-2.127 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


EF/OF 


QO. 

0.337 
0.317 
0.324 
0.306 
0.337 
0.375 
0.259 
0.279 
0.271 
0.386 
0.379 


ll 
0.386 
0.259 
0.126 
0.325 
0.044 
0.013 

13.457 


0.621 


0.529 
0.125 


-1.555 
9.537 


22473 
7.618 


11.663 


41.710 
-12.095 


(REDUCED 


HI/AG 


0.529 
0.552 
0.557 
0.516 
0.549 
0.569 
0.565 
0.569 
0.538 
0.532 
0.507 
0.558 


12 
0.569 
0.507 
0.063 
0.545 
0.021 
0.906 
3-811 


MAJOR 


0.937 


0.573 
0.058 


0.694 
32685 


132442 
24-667 


48.018 


-1.943 
-0.680 


BLVARIATE ANALYSES (REDUCED MAJOR 


A 
SIG(A) 


B 
D(RMA) 


MEAN Y 
MEAN XK 


VMAX J 


top 1 
BASE IL 


0.701 


0.285 
0.059 


-0.847 
7.839 


5.967 
23.925 


20.512 


2-542 
-2-127 


9-621 


0.529 
0.125 


-1.555 
9.537 


2-473 
7.618 


11.663 


41.710 
-12.095 


0.937 


0.573 
0.058 


0.694% 
3.685 


13.442 
24.667 


48.018 


-1.943 
-0.680 


AND LAKE WORTH, FL 


LosaM 


0.327 
0.327 
0.329 
0.315 
0.331 
0.379 
0.299 
0.3396 
0.294 
0.304 
0.299 
0.293 


12 
0.396 
0.293 
0.103 
0.324 
0.033 
0.009 

10.127 


AXIS) OF Y/X. 


0.807 


0.514 
0.088 


4.597 
6.673 


7.783 
23.925 


34-020 
15.115 
-B 2494 


AXIS) OF LOG 


0.835 


1.614 
0.256 


-1.336 
2.929 


0.886 
1.377 


1.690 


5.584 
-27.116 


Mey MCZ 1979424232557,232553 -- LIVING 


LO/AG 


0.312 
0.319 
0.326 
0.299 
0.326 
0.358 
0.302 
0.377 
0.296 
0.287 
0.289 
0.284 


12 
C.377T 
0.284 
0.093 
0.314 
0.029 
0.008 
9.189 


0.055 
0.053 
0.053 
0.055 
0.053 
0.055 
0.055 
0.065 
0.053 
0.055 
0.048 
0.060 


12 
0.065 
0.048 
0.016 
0.055 
0.004 
0.001 
72185 


RIBHT/AG RIBHT/AP 


0.057 
0.055 
0.054 
0.058 
0.054 
0.059 
0.054 
U.068 
0.053 
0.059 
0.2050 
0.962 


12 
0.068 
0.050 
c.018 
0.057 
0.005 
u.-00L 
8.204 


RIB/GRV 


1.514 
1.519 
1.536 
1.400 
1.300 
2.368 
1.200 
1.905 
1.750 
1.667 
1.458 
1.611 


12 
2.368 
1.200 
1.168 
1.602 
0.307 
0.089 

19.178 


RIB/AM 


0.075 
0.056 
0.062 
0.065 
0.061 
0.070 
0.047 
0.063 
0.058 
0.059 
0.060 
0.056 


12 
0.075 
0.047 
0.028 
0.061 
9.007 
0.002 

11.702 


RESULTS NOT VALID FOK COLUMNS OF LOG Y/LOGX. 


0.850 


0.487 
0.074 


4.235 
5.961 


7.783 
24.667 


37.179 


21.068 
-7.981 


Y/LOG X, VALID ONLY 


0.850 


0.487 
0.074 


4.235 
5.961 


7.783 
24.667 


37.179 


21.068 
-7.981 


0.806 


0.062 
0.011 


-0.165 
6.393 


1.358 
24.667 


5.085 


22.191 
-2-314 


0.806 


0.062 
0.011 


-0.165 
6.393 


1.358 
24.667 


5.085 


22.191 
-2.314 


0.741 


0.065 
0.013 


-G.200 
7.230 


1.358 
23.925 


4.685 


10.585 
-2.797 


0.675 


1.291 
0.275 


0.257 
15.990 


1.464% 
0.935 


6.714 


10.817 
2-292 


FOR COLUMNS 5 


0.741 


0.065 
0.013 


-0.200 
7.2230 


1.358 
234925 


4.685 


16.585 
-2.797 


0.675 


1.291 
0.275 


0.257 
15.990 


1.464 
0.935 


6.714 


10.817 
2.292 


0.790 


O.113 
0.020 


-1.248 
6.543 


1.464 
23.925 


7.255 


15.632 
-11.120 


0.790 


O.113 
0.020 


-1.248 
6.543 


1.464 
23.925 


7.255 


15.632 
~11.120 


(A0)/ Am 


0.055 
0.055 
0.055 
0.054 
0.053 
0.058 
0.052 
0.065 
0.048 
0.044 
0.048 
0.052 


12 
0.065 
0.044 
0.021 
0.053 
0.005 
0.002 

10.081 


0.828 


0.086 
0.014 


-0.781 
5.906 


1-274 
23.925 


5.660 


1.429 
6.954 


0.828 


0.086 
0.014 


-O.781 
5.906 


1.274 
23.925 


5.660 


1.429 
-6.954 


(CEI/AG 


0.059 
0.055 
0.063 
0.060 
0.060 
0.077 
0.060 
0.069 
0.052 
0.051 
0.053 
0.053 


12 
0.077 
0.051 
0.026 
22059 
2.008 
2.002 

13.138 


0.743 


0.108 
0.021 


-1.168 
7.397 


1.464 
24-667 


7.952 


10.943 
8.957 


0.743 


0.108 
9.021 


-1.188 
72397 


1.464 
24-667 


7.952 


10.943 
—B.957 


35 


36 


AEQUIPECTEN NUCLEUS 
12 RIGHT VALVES 


SPEC (AD)/ (CE) 


le 
4. 
66 
qT. 
26 
ll. 
De 
12. 
10. 
8. 
3. 
9. 


NUM. 
MAX. 
MIN. 
RANGE 
MEAN 
STDDV 
SDOVMN 
v 


0.891 
0.976 
0.86% 
0.854 
0.872 
0.712 
0.868 
0.891 
0.935 
0.812 
0.875 
0.964% 


12 
0-976 
0.712 
0.26% 
0.876 
0.079 
0.02) 
7.962 


(IK)/AG 


-O. 
0. 
-0. 
“0. 
-0. 
0.336 
-O6~ 
-0O. 
0.293 
-0. 
-O. 
O06. 


2 
0.336 
3.293 
0.043 
0.315 
0.030 
0.021 
9.581 
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(PQ)/AM (IK)/(PQ) 


-0. 
-0. 
-0O. 
-O. 
-O. 
0.276 
-0. 
-0. 
0.239 
-0. 
-O. 
-0. 


2 
0.276 
0.239 
0.037 
0.258 
0.026 
0.019 

10.217 


UIVARITATE ANALYSIS (REDUCED MAJOR 


R 


A 
SIG(A) 


6 
D(RMA) 


MEAN Y 
MEAN X 


YMAX 1 


TOP 1 
BASE 1 


0.875 


0.758 
0.106 


9.165 
8.803 


1.274 
1.46% 


5.989 


-3.663 
1.470 


1.000 


0.679 
0.000 


-8.668 
0.002 


72550 
23.909 


49.011 


392945 
150731 


1.000 


0.923 
0.000 


-15.462 
0.001 


6.000 
23-2250 


53.769 


166462 
-36.077 


O. 
O. 
0. 
O. 
OQ. 
l. 
0. 
0. 
le 
0. 
0. 
0. 


l. 
le 
Oo. 
1. 
O. 
0. 
36 


AXIS) OF Y/X. 


le. 


le 
0. 


=e 


0. 


Te 
66 


31. 


Te 
-3. 


-- MIAMI AND LAKE WORTH, 


288 


222 


2 
288 
222 
066 
255 
046 
033 
699 


009 


583 
000 


950 
093 


550 
000 


300 


804 
539 


FLAse, MCZ 197942,232557,232553 -- LIVING 
(BJ)/AM (AB)/AG (BJ)/(AB) (MN) /AM 
-0. -06. O. -0. 
-O. -0. QO. -O. 
-O. -0. O- -O. 
-0. -O~ O-. -0. 
-0O. -0. Ue -0O. 

0.397 0.194 1.939 0.276 
-U. -0. Ue -0. 
-0. -O. O. 0.280 

0.354 0.218 1.633 0.217 
-O. -O. O. -O~ 
-0. -0. Oe 0. 
-O. -O. O. 0.238 

2 2 2 4 

0.397 0.218 1.939 0.280 

0.354 0.194 1.633 Q.217 

0.044 0.024 U2 306 0.063 

0.376 0.206 1.786 0.253 

0.031 0.017 Je2l6 0.030 

0.022 0.012 Oo153 0.015 

8.188 8.284 12.122 12.024 


1.000 


1.154 
0.000 


-18.077 
0.001 


8.750 
23.250 


68.462 


59.213 
-28.573 


QO. 


0. 
0. 


4.909 
11.477 


4.900 
23.900 


4.900 


-32.993 
16.007 


Ue 


QO. 
C. 


8.750 
14.957 


6.750 
4.900 


8.750 


-35.169 
13.831 


0.835 


0.510 
0.140 


=5'. 724 
5-900 


52675 
22.350 


322528 


30.693 
-14.024 


(FG)/AG 


-0O. 
-O. 
“VU. 
0.094 
-0. 
-0. 
-0. 
0.109 
0.076 
-0. 
0.084 
0.106 


5 
0.109 
UV-0T6 
0-033 
0.093 
0.014 
0.006 

lo.1lly 


RESULTS NOT VALID FOR COLUMNS OF LOG Y/LOGX. 


0.589 


0.178 
0.064 


-1.957 


8.850 


2.160 
23.100 


13.191 


22.819 
-10.653 


AK/GP 


1.154 
1.071 
1.110 
1.103 
1.164 
1.033 
1.034 
1.171 
1.198 
1.095 
1.119 
1.071 


12 
LeL7l 
1.033 
0.138 
1.192 
0.045 
0.013 
4-122 


9.921 


1.146 
0.129 


-0.516 


4.146 


12.817 
11.633 


45.328 


6.527 
-0.633 


AK/AL) 


1.014 
1.063 
1.044 
0.97) 
1.024 
32992 
329927 
1.066 
L.07)9 
JieD2il 
1.034 
1.05u 


12 
1.07, 
3-927 
vel 42 
1.02> 
2004S 
2.013 
42264 


Je9L7 


1.069 
Be l'22 


=)2961 
42224 


17-817 
12.567 


412344 


8.99T 
-).68H 


AEQUIPECTEN NUCLEUS 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 


12 RIGHT VALVES 


SPEC €F/0G 
Le Oo. 
4. 0.216 
6. 0.205 
Te 0.195 
2. 0.185 
Ale 0.228 
Se 0.259 
12. 0.167 
10. 0.173 
8. O.171 
3e 22245 
9- e253 
NUM. ll 
MAX. 0.259 
MIN. 0.167 
RANGE 0.092 
MEAN 0.209 
sToov 0.034% 
SOVAN 0.019 
v 16.202 


LOG PA/AM 


0.083 
0.041 
0.048 
0.001 
0.024 
0.093 
-0.042 
0.087 
-0.149 
-0.175 
-0.152 
-0.216 


12 

0.093 
-0.216 
0.309 
-0.034 
0.112 
0.032 
—333.033 


PA/AM 


0.047 
0.042 
0.045 
0.041 
0.039 
0.056 
0.037 
0.056 
0.928 
0.026 
0.029 
0.027 


12 
0.056 
0.026 
0.030 
0.039 
0.011 
0.003 

26.883 


BIVARIATE ANALYSIS (REDUCED MAJOR 


R 0.388 
A 0.370 
SIG(A) 0.103 
8 -1.925 


DIRMA) 10.767 


MEAN Y 2.473 
MEAN X 11.882 


YMAX 1 12.880 


TOP 1 51.179 
BASE 1-14.973 


0.802 


34426 
0.590 


~4.769 
7el71 


-0.043 
1.377 


1.664 


13.689 
-179.271 


0.761 


0.128 
0.024 


-2.112 
6.990 


0.955 
23.925 


72503 


-30.618 
-S.174% 


BIVARIATE ANALYSIS (REGRESSION OF 


R 0.388 
AYX) 0.144 
BtY) 0.765 
A(Kky) 1.049 
Bix) 9.287 
St¥.x) 1.021 
S(K.Y) 0.478 


D(y.X) 15.283 


S(DYX) 3.258 
D(X.Y) 8.593 
S(DXY) 1.832 
CONAYX 9.114 
CONBY 9-129 
CONAXY 0.830 
CONBX 0.325 
YMAX 2 6.514 
AMAX 15.897 
TOP YX 1.667 
BASEYX 5.949 


TOP XY¥-21.099 
BASEXY 8.127 


0.802 


2.749 
-3.828 


0.234 
1.387 


0.026 
0.090 


-210.843 
-43.038 
1.908 
0.390 


0.647 
0.027 
0.055 
0.008 


1.327 
1.692 


0.982 
-144.169 
-3-427 
25.886 


0.761 


0.098 
-1.379 


52937 
18.253 


1.558 
0.200 


20.905 
4-267 
6.512 
1.329 


0.026 
0.060 
1.600 
0.472 


5.939 
137.000 


-34.450 
-3.379 
28.933 

8.518 


~~ MIAMI AND LAKE WORTH, FL 


OF/AM 


-0. 
0.325 
0.315 
0.305 
0.300 
0.360 
0.335 
0.343 
0.301 
0.315 
0.320 
0.342 


LL 
0.360 
0.300 
0.060 
0.32% 
0.019 
0.006 
5-982 


AXIS) OF Y/XK. 


0.815 


0.364 
0.064 


-0.961 
5.562 


7.618 
23.555 


26.357 


2-659 
-1.884 


yY ON XK AND K 


0.815 


0.297 
0.622 


2-239 
6.499 


1.233 
0.449 


5.893 
1.256 
52233 
1.116 


0.070 
0.143 
0.530 
0.392 


22.899 
62.468 


-G.119 
1.219 
-5.848 
3.033 


Aes MCZ 197942,232557,5232553 -- LIVING 


LO/ALM 


0.269 
0.267 
0.271 
0.264 
0.274 
0.283 
0.251 
0.292 
0.246 
0.250 
0.252 
0.257 


12 
0.292 
0.246 
0.046 
0.265 
0.014 
0.004 
5.361 


RESULTS NOT VALID FOR COLUMNS OF LOG Y/LOGX. 


0.964 


0.353 
0.027 


2.544 
3-292 


7.783 
29.292 


32.712 


12.650 
4.7194 


ON Y) OF Y/X. 


0.964% 


0.340 
-2-170 


2.734 
8.013 


0.932 
0.329 


4.223 
0.862 
3.183 
0.650 


0.030 
0.100 
0.239 
0.282 


31.811 
79.094 


10.952 
-4.091 
Te Slit 
2.805 


OG/GP 


1.115 
0.986 
1.000 
1.051 
1.082 
1.050 
1.000 
1.135 
0.991 
1.057 
1.009 
1.010 


12 
1.135 
0.986 
0.149 
1.041 
0.050 
0.014 
4.7194 


0.891 


1.095 
0.144 


-0.644 
4.900 


12.100 
11.633 


43.174 


3.688 
-0.789 


0.891 


0.976 
0.750 


0.813 
1.796 


0.539 
0.591 


4.881 
0.996 
4.635 
0.946 


0.158 
0.179 
0.131 
0.163 


39.776 
34-317 


-0.275 
0.919 
-4.973 
1.571 


P.AR/LO 


0.144% 
0.128 
0.137 
0.131 
O.117 
0.149 
0.123 
O.141 
0.094 
0.086 
0.096 
0.093 


12 
0.149 
0.086 
0.062 
0.120 
0.022 
0.006 

16.509 


0.970 


0.250 
O.OL7 


-0.987 


32949 


0.955 
7.783 


5.501 


-35.522 
2.418 


RESULTS NOT 


0.970 


0.242 
-0.929 


3.888 
4.069 


0.298 
0.074 


7.791 
1.590 
3.832 
0.782 


0.019 
0.023 
0.306 
0.090 


32366 
81.836 


-35.852 
-2.277 
752164 

5.417 


OF/ALM 


-O. 
0.266 
0.260 
0.255 
0.248 
0.269 
0.2861 
0.253 
0.252 
0.259 
0.269 
0.300 


1 
0.300 
0.248 
0.052 
0.265 
0.015 
0.005 
5.663 


0.903 


0.233 
0.030 


0.885 
5.048 


7.618 
28.864 


24.214 


“1.541 
1.734 


VALID FOR 


0.903 


0.211 
1.540 


3.869 
-O.611 


1.431 
0.334 


4.383 
0.934 
4.959 
1.057 


0.033 
0.106 
0.615 
0.455 


22.597 
96.114 


-4.710 

3.019 
-1.360 
-0.214 


0. 


0. 
oO. 


0. 


0. 


oO. 
oO. 


oO. 


oO. 
oO. 


COLUMNS OF LOG Y/LOG X. 


0. 


QO. 
oO. 


QO. 
Oo. 


oO. 
0. 


QO. 
QO. 


QO. 
QO. 
QO. 
QO. 


9. 


oO. 
0. 


0. 
Oo. 


0. 
0. 


0. 
oO. 
oO. 
QO. 


0. 
oO. 
Oo. 
0. 


0. 
oO. 


QO. 
oO. 
0. 
oO. 
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Error Messages in DASAN 


Only one error message has been written into DASAN, 
which reads, ANOMALOUS VALUE IN DATA. This 
statement will appear and the job will be terminated for 
any of the following reasons: 

1. Control cards contain negative values where these 
are not allowed. 


2. A zero on a card in the data deck has been used as 
a signal for division by a constant, whereas only a +1 ora 
—| is permitted. 

3. A variable which is to be converted to a logarithm 
is negative. 


VPLOT: Description of Program 


The variable plotting program (VPLOT) provides for very 
rapid plotting of bivariate scatter diagrams and the con- 
struction of reduced major axes or regression lines. The 
univariate output array of DASAN (Figure 1), or any array 
in which columns represent variables (measurements) and 
rows represent specimens, may serve as the input data array. 
Within VPLOT, this input array is altered and/or expanded 
to form an internal data array. Alterations consist of loga- 
rithmic (base 10) transformations of any of the input vari- 
ables, with the logs replacing the original values, Alterna- 
tively, the logs may be entered in new columns of variables, 
thereby preserving the original values. An alteration of the 
input array is desirable if only the logs are to be plotted and 
not the original values. On the other hand, expanding the 
input array by entering the logs in new columns allows both 
the original data and the logs to be plotted. 

The main routine of VWPLOT is concerned with input, 
with the alteration and expansion of the input array, and 
with output. The actual plotting of each bivariate scatter 
diagram is carried out by the subroutine, APLOT, which 
employs both Boolean algebra and FAP. This subroutine 
accepts and plots the values in two columns of the internal 
data array at a time according to column indentification 
numbers specified on control cards. 

Listed below are the main steps followed by the program: 


VPLOT (main) 


1. Control cards, variable format, and data are read in. 

2. The value of each input variable (except for the first, 
which should be a specimen indentification number) is 
compared to a specified minimum value. If below this value, 
the variable is replaced by the specified minimum. 


3. The value of each input variable is compared to a 
specified maximum value. If greater than this value, the 
variable is replaced by the specified maximum. 

4. Columns of variables, specified on control cards, are 
replaced by their logarithms (base 10). 

5. The logarithms of specified columns of variables 
are computed and entered as new columns of data. 

6. Two variables, designated as x and y on control cards, 
are entered as arguments in a call to Subroutine APLOT. 


Subroutine APLOT 


7. One bivariate scatter diagrams is computed, and 
control is returned to the main program. 


VPLOT (main) 


8. The scatter diagram computed by APLOT is 
printed. 

9. If more diagrams remain to be plotted, Steps 6 
through 8 are repeated for each plot. 

10. If more data decks remain in the run, control returns 
to Step 1 with the reading in of new control cards and data. 
Otherwise, the job is completed. 

Computation Time.—The sample problem described be- 
low (12 specimens, 36 input variables, and 43 plots) required 
only 50 seconds on the IBM 7094. Increasing the number of 
specimens does not add greatly to the computation time. 

Core Space.—The space occupied by VPLOT is not large, 
but because there are difficulties in plotting points consisting 
of three digits rather than two, 99 must remain the upper 
limit of the number of specimens which can be run at a time. 
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Listing of FORTRAN II Statements in VPLOT 
VARIABLE PLOTTING PROGRAM 


THOMAS R. WALLER, DEPT. OF PALEOBIOLOGY, SMITHSONIAN INSTITUTION 


NOW 


4d 


HS 


14 


15 
16 
43 


FEBRUARY 6, 1967 


DIMENSION X(99,100),A(99),B(99),TAG(100), IDYVAR(20) , IDXVAR(10),XA(VMAIN 1 
15), YA(5),FME (60), ANAME(15 ) ,LALCOL (60) , LENCOL (60) 


READ INPUT TAPE 5,37,NVAR,NTAG 
READ INPUT TAPE 5,32, (TAG(J),J=1,NTAG) 


READ INPUT TAPE 5, 32,ANAME 

IF (ANAME(1)-(+6HFINISH) )2,40,2 
Tet 530 

READ INPUT TAPE 5,37,NS 

MIN = NS +1 

NS2 = NS + 2 

NS3 = NS + 3 

K = 
READ INPUT TAPE 5,37,NGROUP 

READ INPUT TAPE 5,31, IDYVAR, IDXVAR,NY,NX,NODATA 
IF (NODATA) 30,4, 24 

READ INPUT TAPE 5,37,NFMI 

KK = 12 * NEM 

READ INPUT TAPE 5, 32, (FME(I),I=1,kK) 

READ INPUT TAPE 5,FMT,((X(I,J),J=1,NVAR), I=1,NS2) 
READ INPUT TAPE 5,37,LALTER,LENTER 

IF (LALTER)30,6,5 

READ INPUT TAPE 5,39, (LALCOL(J) ,J=1,LALTER) 

IF (LENTER)30,41,7 

READ INPUT TAPE 5,39, (LENCOL(J) ,J=1,LENTER) 

NVARL = NVAR + 1 

NVAR2 = NVAR + LENTER 

IF (LALTER + LENTER)30,8,45 

READ INPUT TAPE 5,42, (X(NS3,J),J=1,NVAR) 


OPN ell 


IF (LALTER)30,18,13 

DO 43 JL=1,LALTER 

J = LALCOL(JL) 

DO 16 I=1,NS2 

IF (X(I,J))30,15,14 
X(I,J) = LOGLOF(X(I,J)) 
GO TO 16 

X(I,J) = X(NS3,J) 
CONTINUE 

CONTINUE 
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18 IF (LENTER)30, 8,19 
19 J2 = NVAR 


20 
el 


22 
23 


10 
46 


aa 
12 
47 


eh 


25 


26 


27 
28 


29 


30 


DO 23 JL=1,LENTER 
= LENCOL(JL) 
Je+l 


)30,21,20 
KOE ae LOGLOF(X(I,J)) 
GO TO 22 
X(I,J2) = X(NS3,J) 
CONTINUE 
CONTINUE 


DO 46 I=1,NS 

DO 10 J=2,NVAR 

IF (X(I,J)-X(MIN,J))9,10,10 
K(I,J) = x(MIN, J) 

CONTINUE 

CONTINUE 


DO 47 I=1,MIN 
DO 12 J=2,NVAR 
mm (x(@,d) — X(NS2,J)))l2,12, 01 


X(I,J) = X(NS2,J) 
CONTINUE 
CONTINUE 

XA(1) = 0.0 
YA91) = 0.0 

i = The 2 

DO 28 M=1,NX 

L = IDXVAR(M) 

DO 27 J=1,NY 


IF (JJ - L)25,27,25 
K=K+1 

DO 26 I=1,NS2 

B(I) = ee 

A(z) = X(I,L) 


WRITE OUTPUT TAPE 6, 34,NS,ANAME 
WRITE OUTPUT TAPE 6,35,K,TAG(JJ) 
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WRITE OUTPUT TAPE 6,36,TAG(JJ),TAG(L) ,TAG(L) 


CALL APLOT (1,NS2,1,1,XA,YA,A,8B) 
CONTINUE 
CONTINUE 


IF (NGROUP - II)29,29,3 
WRITE OUTPUT TAPE 6, 33 
GO TO 1 

WRITE OUTPUT TAPE 6,38 


DuoFrwrmr 
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32 FORMAT (12A6) 

33 FORMAT (1H1/1H1) 

34 FORMAT (1H1,10H PLOT FOR ,1T4+,16H SPECIMENS » 15A6) 
35 FORMAT Ree NO. , I3,40X,11HORDINATE - ,A6) 
36 FORMAT (1HO,5X,A6,8H VERSUS ,A6,32X,11HABSCISSA --,A6) 
37 FORMAT (I3) 

38 FORMAT (1H1,24HNEGATIVE VALUE IN INPUT.) 

39 FORMAT (3612) 

42 FORMAT (F4.0,9F7.4/(1LOF7.4) ) 


31 FORMAT fey 


40 CALL EXIT 
END 


SUBROUTINE APLOT(LeNaNXAXISSNKYAXIS sXAdYASX9Y) 


DIMENSTON X(140) 9¥(140) »sPLOT (50912) »XA(5) sYA(5) s0UT(&8) 95 (140) 9ALPn 
1A(275) 


FORMAT(1H »84X5HALPHAs10X14HIDENTIFICATION) 
FORMAT(1H s6XlHesl2A6s9lHeslll(2XxA2) ) 

FORMAT (1H 9TXaT(lLOH+Veccecccee) sliitsliieslHt) 
FORMAT (1H 94XF60294XF60294XF 602 94XF 6029 4X6 02 94XF 602 94XF 602 94XF E02 ) 
FORMAT(1H sF6e291H+e912A6591rit+911(2XA2)) 
FORMAT ( 1HOs2HDX 92XF 8e492X2HDY 9 2XF B04) 
LSAV =L 

NSAV = N 

= Xa(GlE,) 

Roms XG) 

We NEC) 

NeSip a Yates) 

Net=se NGG a) 


DO 14I1=L95N 

IF (XL-X(1))75898 

Xs =. XC1) 
LEOXSS XCD) 10's 10's: 2 
xSie= xl) 

MiG VAI) Lalas 12h) 2 
VERS yale) 
LECYS=V¥ CI) Waal 45 13 
¥Se= Ye) 

CONTINUE 


ITF (XE—-XS)1029102915 
HE CVE—V¥ Sy VO 2 sO2 16 
DX (XL=XS)/35e0 
DY GY E=Y¥SV749'e.0 


DOWZOONK = ts5:0 
DO 300 J WoRleZ 


42 


300 PLOT(KsJ) 
DOP 200 1 

ALPHA(T) 

ALPHA(1) 

ALPHA(12) 
ALPHA( 23) 
ALPHA (34) 
ALPHA(45) 
ALPHA(56) 
ALPHA(67) 
ALPHA( 78) 
ALPHA( 89) 
ALPHA(100) 
ALPHA(111) 
ALPHA(122) 
ALPHA( 133) 
ALPHA(144) 
ALPHA(155) 
ALPHA(166) 
ALPHA(177) 
ALPHA(188) 
ALPHA(199) 
ALPHA(210) 
ALPHA(221) 


200 


HT eit 
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(+6H ) 


L275 
(+6H ) 


(+6HAAAAAA) 
(+6HBEBEbB) 
(+6HCCCCCC) 
(+6HUDDDDV) 
(+6HEEEEEE) 
(+6HFFFFFF) 
( +6HGGGGGO) 
(+6 HAHHHHH ) 
(+6H9 Jd) 
(+6HKKKKKK) 
(+6HLLLLLL) 
(+6HMMMMMi4 ) 
(+6HNNNNNN) 
(+6HOOO0OO) 
(+6HPPPPPP) 
(+6HOGQGQQ) 
(+6HRRRRRR) 
(+6HSSSSSS) 
(+6HTTTTTT) 
( +6HUUUUUU) 
(+6HVVVVVV) 


(+6HwWWwwWWw ) 
(+6HXXXXXX ) 
ALPHA( 254) CEG YY Yo) 
ALPHA( 265) (+6HZ222222Z) 
B(l) 00¥109010001 
BiG) 008200020002 
B(3) 009300030003 
B(4) 006400040004 
Bi!) 0005000500905 
B(6) 000600060006 
b(7) 000700070007 
6(8) C01000100010 
B(9) 00110C110011 
B(10) =010001000100 
Belly =OLOLOLOTOT OL 
B12) =OLC 201020102 
Bilis) =OLRCZOL03 01038 
B(14) =010401940104 
B(15) =015501050105 
B(16) =010601060106 
B(17) =C19701070107 
B(18) =011001100110 
BICHON, =O JO ONe nT 
B(20) =020002000209 
B(21) =020102010201 
B(22) =024202020202 
B(23) =020302030203 
B(24) =02C402040204 


ALPHA (232) 
ALPHA (243) 


CWoODBDWBWOetWowTwWecewTOeocCvwwwWoeococoawwo 


CoOMDWOWBIWBWIAOWOCWOACCAOeUBTOCGoOOCctCuvoOWWOOvewtwvwoowawecwtcnwnaowenwnwwecw wD 
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B(25) =020502050205 
B(26) =020602060629 


BIC2e) sO ZOO COMO ZO 
B(28) =021902100210 

SUZ) =O21TO2TT Oa 1 I 

BBG) T= OZVCCA0COS 00 
B(31) = 030103010301 
B(32) = 030203020302 
Bit33) = 030303030303 
B(34) = 030403040304 
BiW3'5)) = 030503050305 
B(36) = 030603060306 
BiG3i)) = OZ 07030703107, 
B(38) = 031003100310 
BiGs SO) =O ST TOS TOs aT 
B(40) = 040004000400 
B(41) = 040104010401 
B(42) = 040204020402 
B(43) = 040304030403 
B44) = 040404040404 
B(45) = 040504056405 
B(46) = C40604060406 
B(47) = 040704079497 
B(48) = 041004100410 
B(49) = 041104110411 
B(50) = 050095000500 
Bist) = ObOL0s010501 
BDZ) = OSD 205020502 
B(53) = 050365030503 
B(54) = 050405040504 
B(55) = 050505050505 
B(56) = 050605060506 
BCs)! = Ods07 OS 0703/07 
B(58) = 051005100510 
BiGso)s = OSIVOS OST 1 
B(60) = 060006000690 
Bites) = O6O10SC 1060 1 
B(62) = 060206020602 
B(63) = 060306030503 
B(64) = 0600406040604 
B(65) = 060506050605 
B(66) = 060606060606 
B(67) = 060706070607 
B(68) = 061006100610 
BiU69)) = O61 VEGtlOG Ty 
B(70) = O700C7TO00790 
Bile )a = OMCakC xO OROn! 
Bi2:) = O79 20102 07.02 
Bigs) = O03 CVOZ0708 
B(74) = 070407040704 
BiGione = Or OSwil Os Oi OS 
B(76) = 070607060706 
Bia? (ae e=) Ol OmOO.7 OO 


298-840 O—69——4 


UNITED STATES NATIONAL MUSEUM BULLETIN 285 


HS 
HS 


B B(78) = 071007100710 
B By TAS) )) SO 7eab alten y/ al ale) 7 al al 
B B(80) = 100010001000 
B B(8ay)y = LOO LNOOLIO OL 
B B82) = LOOZVOO2T 002 
B B(83) = 100310031003 
B B(84) = 100410041004 
B B(85) = 160510051005 
B B(86) = 100610061C06 
B B(87) = 100710071007 
B B(88) = 101010101910 
5 Bri) s=F UGA LOMO 
B B(90) =110011001100 

B BICSHS) =s NO OMmle OM: 
3 Bag 2 Fr = NO TOA O2 
B B93) = LOS LLOS TL NO3 
B B(94) = 110411041104 
B BICC Sy) = s LOSI OS hos 
B B(96) = 110611061106 
B BICS PSP L oT UC rE Oy 

B BIS) = DAE WOL TNO PRO 
3 Bi y= ea eee ede 


IF (NXAXIS) 32932917 


G 
17 DO 31 1 = l»eNXAXIS 
IF (XA(1)-XS)31918918 
18 IF(XA(1)-XL)19919931 
19 DIFX = XA(L)-XS 
DEUX== Dx 
el 
ZO MIROS! = 91 
21 IF (DIFX-DELX) 24922522 
22 DELX = DELX+DX 
IPOS = IPOS+2 
IF (1PO0S-5)21921923 
23 J = Jtl 
GO TO 20 
24 GO TO (25925927927929929) s1POS 
DRDOn 26 =) si50 
26 PLOT(Ksd) = (46H ) 
GO TO 31 
PUI 28 Ke — sb 0 
ZO PEOT UKs), = (+6H. 91) 
GO TO 31 
29EDOI30°K = 1s 50 
30 PLOT(KsJ) = (+6H 1) 
31 CONTINUE 
€ 


32 IF(NYAXIS) 41941933 

33 DO 40 I = IsNYAXIS 
IF(YA(1T)-Y¥S) 40934934 

34 IFCYA(1I)-YL) 35935940 


DWDOOOUODODW @ 


@ 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 45 


35) [DAN YC SONIC CON ESS 
DEE == DY 
Keel 


Sor UR CONE YS DiS) 389379 37 
STD Ela Ye = DEY DY, 

K = K+l 

CKO IO 2x5) 
38 DO 3 9) =Aserl2 
BSE PROTEC) = eh ==———— ) 
40 CONTINUE 


LADO, SS. = ieisiN 


DiElLDG -yD).4 

DBE a=" DY, 

Dix t= Xe (le) = XcS) 
Dee = alias 
J= 1 

KS 


42 IF(DIFY-DcLY)44943 943 
G3) DE EW] = IDiEBEYA DY. 

Ky = Kat 

GOTO 42 
ae BP OS: ss al 
45 IF (DIFX-DELX) 48 946946 
46 DELX = DELX+DX 

IPOS = ShPROS+2 

IF (IT PO0S-5)45345547 


Qf eo = 
GO TO 44 

48 GO TO (49349362962575375)sI1POS 

49 D = (6060000000C0+(77770000CCNO*PLOT(KsJ)))*(-—(606000000000*PLOT(K 
Ney) 
E = (404000C90000+(777729000000C*PLOT(KsJ)))*(-(404000000000*PLOT(K 
9iJ)))) 
F = (003100000000+(0077000CQ000*PLOT(KoJ)))*(-(0031000000900*PLOT(K 
Le EY 


G = 60000C00COXPLOT(KsJ) 
Li GD)F510%s)5:3 i550 

IOP SE GED as 5'3 5 5:1 

Se GR DSi) oiZ 

52u sR Gy) Siiseo1 4915) 7 


SS PEOTCKs) = (OOOO T li it tse Ol Ch s J) MEG Tai HO COOOOOC#BiCI)>) 
GO TO 88 
54 DO 55 M = 239266911 
IF (ALPHA (M)—(+6H ) )515's\ 516151515 
55) ‘CONT ENUE 
GO 1:04 '88 
56 ALPHA(M) = PLOT(K»sJ) 
ALPHA(M+1) = B(1) 
PEON) =(COOCTIIT17 7 C#HPLOT (Ke J) C77 7TTOOOOCOOOXALPHAIM=1))) 
GOP LON 88 
27 DO! S56) Mi = Ts265%5 11 
IF ( (77770CQOO0000*PLOT (Kod) )*(777700000000% (—(ALPHA(M)))))58959958 


46 


@D 


JB 


OwowodcdD 
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58 CONTINUE 
GO TO 88 
59 P= MZ 
KK = M+10 
DO 60 MM = JJsKK 
IF (ALPHA (iM) —-(+6H ))609615s60 
60 CONTINUE 
GO TO 88 
61 ALPHA(MM) = BCI) 
GO TO 88 
62 D = (O0006C690CO00+(00COT7TTTOOOC#PLOT (KJ) ))*(-(O00006E0600000*PLCT(K 
IED) 
E = (0000404000090+(00N07T7T7TTOONOO*PLOT(KsJ)))*(-(O00040400000*PLCT(K 
Les) 
F = (0000003190004 (0C20909077TIONOXPLOT (K9J)))*(-( 0000003 10000*PLOTI(K 
I) 
G = OOOCECOCCCODO*PLOT(KsJ) 
IF (D)63966963 
63 TEWEN64s65 564 
64 IF(F)65966365 
65> TCG) 70is:6 7% 57/0 


66 REO CKs Sy) = rita OOOO Hl Te REOT EK 63) HU D0O0AT 11-0000 Xb CL)) 
GO TO 88 
67 DO 68 M = 29266911 
IF (ALPHA(*\) —(+6H )16385s69968 
68 CONTINUE 
GO TO 88 
69 A = PLOT(KsJ) 
ER 
CALL MOVE(AsE91251) 
ALPHA(™) = E 
ALPHA(M+1)= 8(T) 
PLOT(KsJ) = (777700007777*PLOT(KsJ))+(000077770000*ALPHA(M=-1) ) 
GO TO 88 


Ol DOM ilar = 9 Ais 266 S)5)-1) 4: 
IFC (OOOS7TTTTOCOO*XPLOT (KsJ))* (000077770000* (-—(ALPHA(M))))) 71972971 
71 CONTINUE 
GO TO 88 
V2 = Mice 2 
KK = M+10 
DO 73 MM = JJsKK 
IF (ALPHA (MM)-(+6H i aeeeoucse rf} 
73 CONTINUE 
GOTO 8s 
74 ALPHA(MM) = B(1) 
GO TO 88 
75 D = (00C0CG0006069+(00CD0D0007TTITTEPLOT(K9J)))*(-(O00000006060*PLOT(K 
Lisi) ) 
E = (0000000040404+(00C0CO0007T7T7T7*¥PLOT(KsJ)))*(-(000000004040*PLOT(K 
KI) D) 
F = (000009000031+(000ODRQ0007TT¥PLOT(KsJ)))*(-—(000000000031*PLOT(K 
o)))) 
G = OOOCOLCOND60D00*PLOT(KsJ) 


76 
a 
78 
a) 
80 
81 


82 


83 


84 


85 


86 


87 


89 
90 
91 


92 
33 


94 
5) 


96 
oni 


TWO FORTRAN 


LGD e7iGiset 9's). 

Fale) lanrvemn oust 

NG Sion 9's1.8 

IF (G) 835980583 
PiGOMe eh) ean ierhararin 
GO TO 88 

DO 81 M = 23266911 
IF (ALPHA(™) —(+6H 
CONTINUE 

GOMT.O) 7sis 

A= POM CK ss) 

B= 1 

CALL MOVE(AsEs24 51 
ALPHA(M) = E 
ALPHA(M+1) = B(TI) 
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T7TTCOOO*#PLOT(KsJ))+(O0000CC9I07T777¥B(1)) 


))81582581 


) 


PLOT(K9J) = (77777777CO00*PLOT(KsJ)) +(000000007777*ALPHA(M-1) ) 


GO TO 88 

DOM Mer= isi Zod 
IF ( (O00C00007777*P 
CONTINUE 

GOT 10) 88 

J ea Mit2. 

KK = M+10 

DO 86 MM = JJsKK 
IF (ALPHA (MM) -(+6H 
CONTINUE 

GOMTO38 

ALPHA(MM) = B(T) 
CONTINUE 


WRT E OUTRUT TAPE 
= 50 

Male 

41 

400 


NOM 


F (M=266) 90390952 
IF (ALPHA (i) -(+6H 
MM =M-1 

NN = M+9 

Wie ROUmEUiie WARE 
GO TO 93 

WREKE OUTPUT MARE 
M = M+11 

Kes a: 
IF(K)10C9100994 
IF (K-L)96995396 
P= VS+ (Baby) 


K 
p 
E 
B 
M 
I 


EEO 
B = B-19ec 
GON TOR 89 


IF UM=266) 9.799 7999 
IF (ALPHA (id) -(+6H 


LOT(KsJ))*(000000007777*(-(ALPHA(M)))))84985984 


))86387986 


VIVG1Ls92991 


6959Psa(PLOT(KsJ) sJ=1l912)9(ALPHA( 1) » T=MMoNN) 


6959Pe(PLOT(KsJ) sJ=1912) 


1)98s99398 


48 


98 


29 


100 


101 


102 


MOVE 


HERE 


NEXT 


AA 
NPL 


RIE 


TEMP 


AR 
NPRT 
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MM Ma 
NN M+9 
WRITE OUTPUT TAPE 6929(PLUT(K9J) sJ=1912) 9 (ALPHA(1) I =MiMaNN) 
GO O93 


WRITE OUTPUT TAPE 692s(PLOT(KosJ) sJ=1912) 
GO TO" 93 

VRE COURE Ui TAPES Gis. 

A = 0e9 

DO 101 M = 198 

OUT(M) = XS+(A*DX) 

A = At+5e0 


WRITE OUTPUT TAPE 694 9(OUT(#) shi=l98) 
WRITE OUTPUT TAPE 696sDXsDY 


Ey =) SAV 

N = NSAV 
RETURN 

END 

FAP 

GOUNT, 310 
ENTRY MOVE 
CLA 294 
STA BB 
CLA 494 
STA HERE 
CLA 

Ze RITE 
CLA 154 
STA AA 
GSA 394 
STA NEXT 
CLA 

ARS 18 
STA NPL 
CAL 

ALS 

TOV BB 
TRA BB 
CLA 154 
STA AR 
CLA 394 
STA TEMP 
GIEA 

ARS 8 
SiA NPRT 
CAL 

ARS 

SLW 

TRA 594 


END 


TWO FORTRAN II PROGRAMS FOR ANALYSIS OF MORPHOMETRIC DATA 49 


Input Data Preparation for VPLOT 


The following section describes the preparation of control 
and data cards for VPLOT, using the same conventions 
as in the previous section describing the input for DASAN. 
Unlike the previous section, however, a third group of 
control cards is used. These Group III controls follow each 
data deck in the run and, like the Group II controls, con- 
tain infromation which affects only the data deck which 
they accompany. The major groups of cards appearing in 
a VPLOT run which contains two data decks are ordered 


as follows: 


Program Source Deck 

Group-I Control Cards 

Group-II Control Cards 

Data Deck 

Group-III Control Cards 

Group-II Control Cards (for following data) 
Data Deck 

Group-III Control Cards 

FINISH Card 


Teil) 
Deshi) 
Be (=o) = 
4. (2) 
3. (1) 


Group I Control Cards 


The number of variables in the input data 
array (from 3 to 100), including the column 
of specimen identification numbers. Punch 
in columns 1-3. 
The number of variables in the internal data 
array (from 3 to 100), including the column 
of specimen identification numbers. Punch in 
columns 1-3. 
Names of variables in the internal data array, 
with each name not to exceed six characters. 
Punch in columns 1-72 on each card, 12 
names per card, with each name centered in 
a field of 6 columns, such that the name in 
the nth field is for the nth column in the 
array. 

Group II Control] Cards 


Any name or sample-identification informa- 
tion to be printed at the top of each output 
page. Punch in columns 1-72 of first card 
and columns 1—18 of second card. Both cards 
must be present even if blank. 

The number of specimens (from 3 to 97) in 
the input data array, punched in columns 1-3. 
The maximum number is 97 rather than the 
dimensioned 99, because the cards containing 
the minima (Instruction 11) and maxima 
(Instruction 12) are read in as data for two 
additional specimens. 


6. 


Ine 


(Oe 


(1-5): 


Ge): 


The number of plotting groups (from 1 to 
999). A plotting group is a series of x and y 
variables such that each x is plotted against 
every y, as explained in Instruction 7. Punch 
in columns 1-3. 

Information for setting up a plotting group 
(Instruction 6), to be punched as follows: 
(a) Columns 1-40: The column-identifica- 
tion numbers of the y-variables (from 1 to 20 
numbers), with each number in a field of 2 
card columns. 

(b) Columns 43-62: The column-identifica- 
tion numbers of the x-variables (from 1 to 10 
numbers), with each number in a field of 2 
card columns. 

(c) Columns 67-68: The number of y-vari- 
ables in the plotting group (the number of 
values punched for Instruction 7a). 

(d) Columns 69-70: The number of x-vari- 
ables in the plotting group (the number of 
values punched for Instruction 7b). 

(e) Column 72: Punch a zero if this is the first 
plotting group preceding the data deck; punch 
1 if this is not the first plotting group and it 
follows the data deck. (See Instruction 18.) 
The number of cards (from 1 to 5) which are 
required to specify the format for reading in 
the data deck. Punch in column 3. 

The format, in FORTRAN II, for the read- 
ing in of data cards. (See Instruction 21 in 
the description of DASAN.) Punch in col- 
umns 1—72 of each of the cards needed. 


Data Deck 


See Instructions 21 and 37 in the description 
of DASAN. The specimen numbers, which 
form the first column of the data array and 
which must be read in, are not used by 
VPLOT. Rather, the specimens are assigned 
consecutive numbers in their read-in order, 
beginning with No. 1, and it is these consecu- 
tive read-in numbers which appear as the 
points in the output scatter diagrams. 


Group III Control Cards 


Cards containing a minimum value for each 
variable in the input array (see also Instruc- 
tion 17). The coordinates of the origin of each 
bivariate scatter diagram are determined by 


50 


ae 


Ds 


Listed below are the control and data cards for a test run 
of VPLOT, the output from which is discussed in the fol- 
lowing section. The input data is the univariate output ar- 


(0-2) : 
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the minima specified for the « and y variables 
in each plot. The difference between the min- 
imum and maximum (Instruction 12) speci- 
fied for each variable determines the scale of 
calibration of each scatter diagram in which 
these variables are plotted. (See explanation 
of sample output from VPLOT.) The cards 
are read in according to the same format as 
that used for the reading in of the data deck 
and must therefore be punched in the same 
way as the data cards, with a specimen num- 
ber being entered as the first variable. 

Cards containing the anticipated maximum 
values for each of the variables in the input 
array, as explained in Instruction 29 in the 
description of DASAN. These values, together 
with the minimum values specified for each 
variable (Instruction 11), determine the scale 
to which the axes of each of the scatter dia- 
grams are calibrated. These cards are read in 
according to the same format as that used for 
the reading in of the data deck and must 
therefore be punched in the same format as 
the data cards, with a specimen number be- 
ing entered as the first variable. 

The number of columns of variables (from 0 
to 60) which are to be replaced by their loga- 
rithms (base 10). Punch in columns 2-3, 
The number of columns of variables (from 0 
to 60) which are to be transformed to loga- 
rithms, with the logs being entered as new col- 
umns of data and the original values being 
preserved. Punch in columns 2-3. 

Column identification numbers of the vari- 
ables which are to be replaced by their loga- 
rithms. Punch in fields of 2 columns in col- 
umns |—72 of the first card and in columns 
1—48 of the second card. 

Column identification numbers of the vari- 
ables which are to be transformed into loga- 
rithms and entered as new columns of data in 
the internal array. Punch in fields of 2 col- 


liye 


18. 


19. 


20. 


(0-7): 


umns in columns 1—72 of the first card and in 
columns 1—48 of the second card. 

Cards containing a minimum log value for 
each variable in the input array which will 
be transformed to a logarithm and for which 
a zero minimum was specified in Instruction 
11. This log minimum will be inserted by the 
program into any blank data cells and, like the 
minima for non-log data described in Instruc- 
tion 11, will determine the coordinates of the 
origin and the calibration of the axes of each 
diagram. 

If the minimum specified for a variable in 
Instruction 11 is greater than zero, the log of 
this value will become the log minimum, and 
it is not necessary to specify any new log mini- 
mum here. The maximum value for each log 
variable is the log of its maximum previously 
specified in Instruction 12. 

Punch an imaginary specimen number in 
columns 1—4 and the log minima in fields of 7 
columns in columns 5-67 of the first card and 
in columns 1—70 of each succeeding card. Con- 
secutive fields represent consecutive columns 
in the input array, so that a minimum for the 
7th variable, for example, must be punched 
in the field which is in columns 40-46. 

The decimal point will be placed by the 
program in front of the last four digits of each 
log minimum. 

If more plotting groups remain, enter one in- 
(Instruction 7) for each 
group. Punch as described in Instruction 7, 
placing a 7 in column 72. 

If more data decks remain to be processed 
having the same numbers of variables in the 
input and internal data arrays and the same 
names for the variables, begin again with In- 
struction 4. 

If no more data decks remain, the last card 
consists of the word FINISH punched in col- 
umns 1-6. 


formation card 


Sample Input Data for VPLOT 


ray punched by the preceding test run of DASAN. On the 
far right side of each control card are numbers which refer 
to instructions in the section on input preparation. 
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- DATA 
036 1 
O41 2 
SPEC DG GP AD AK AM DF CD EI BJ DE LO 3 
RIBHT HI PLICAEANTCOSPOSCOSINTRIBINTGRV (AD) (CE) (IK) (PQ) (BJ) 
(AB) (MN) (FG) ALM P.AREA KM EF AG CF LOGPARLOG LOLOG AM 
LOG AMLOG DFLOGPARLOG LOLOGALM 
AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLA., MCZ 197942, 232557, 4 
232553 -- LIVING 
012 5 
016 6 
0607083314121321222527 32 1101 © ? 
1 8 
GAP Ort6. 1 lenS. LOK, MHA /I4ES. 159K, DHB/5SF5.a, 57.4. 55m, LHC) 9 
ea am = = INSERT DATASDECK MERE =\S === = —=9=.=- ==> 5 == =| Ho 
98. a 
11 
11 
Gone OO: +On0445 0240.0) 7550 25.20 2720), 1420) 19, 5525.0 26.0) 545) 5050 12 
FORORMOMOnOUO) 55) Seb! food) 6 > 27. One Is.O. 50. On 15012020). 9. O100).0 
PONOmS5 a0" 6.5) 65.0 55.0 1.5010 1.4150 1675l 
00 1 
05 14 
0607291228 16 
98 -10000 -10000 aly 
-10000 eZ 
-10000 -10000 ay) 
ez 
0910121318202324262907 06 sO aE Tt 18 
040731 02 0301 1 
060712 28 0301 1 
05 30 OnOia 
0805 o4 0201 1 
07 08 0101 1 
Bil 07 OnOaNaD 
0502 03 0201 1 
18 19 Odoly aL 
20 2a ple nies! 
22 23 OnOay a 
2h 25 O1On 2 
39 37 0101 1 
29 12 Oot 
3840 Bi, 0201 1 
FINISH 20 
Sample Output from VPLOT 
Reproduced below are two of the plots resulting from an variables in the input array, 41 variables in the internal 


execution of VPLOT involving 12 specimens, with 36 array, and a total of 43 ratios to be plotted. 
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PLOT FOR 12 SPECIMENS AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLlAce 
PLOT NOw- 22 ORDINATE - OF 
DF VERSUS AM ABSCISSA -- AM 
25.00+ 14+ 


met 





20-41+ ‘ 

15.31+ ‘ 

10.20+ + 

5.10+ + 

0. +1376--f-------- 29 ----- 01 ------- -- ------ $n nr rrr : 
Le fal ote atelele Pelate wialeisieledic e/c\s(cleislaje treieleisials aleleit sje lelelc,eie'e\e\t a/a, 10 c1s)a/e'e t sjelcleieioieies tele: 
0 10.71 21.43 32214 42.86 53.57 64.29 75.00 


DX 261429 DY 0.5102 


285 


MCZ 197942, 


AA 
BB 
cc 


ALPHA 
02 
03 
05 


oT 
08 
12 


232557 232553 -— LIVING 


10 


IDENTIFICATION 
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PLOT FOR 12 SPECIMENS AEQUIPECTEN NUCLEUS -- MIAMI AND LAKE WORTH, FLAey MCZ 1979425 232557, 232553 -- LIVING 
x 
PLOT NO. 43 ORDINATE - LOG LC ¥ 
’ 
LOG LO VERSUS LOG AM ABSCISSA --LOG AW 
ALPHA IDENTIFICATION 
1.41+ I 14+ AA O1 04 06 
; I - BB O05 08 
5 I =) CG ..03\ 10 
: I 5 
S I a 
. I B 
x I a 
A I x 
. I s 
0.97+ I + 
5 I ' 
; I : 
: I 5 
: I 3 
a I e 
3 I . 
5 I 5 
a I : 
: I : 
0.48+ I + 
‘ I e 
J I é 
6 I js 
3 I : 
; I 3 
z I . 
s I - 
5 I a 
; I s 
aQo Oller e Se ose Sama SSS SS Se SS SS SS SSS + ar iain aimiaiaaaiammmeacl + 
; I 5 
‘ I . 
A I 5 
° I 5 
: I : 
: I e 
5 I a 
: I 4 
* I Hi ° 
-0.51+ I Us ca 
. I Vf . 
. I he . 
g I yh ss 
; I . 
; Teme ; 
; [ey . 
; I : 
-1.00+13 I A + 
Pocccccccceteccccccectecccccceetesccccccetocscccccctecccccccctecscccccetet 
-1.00 -0.59 -0.18 /0.23 0.64 1.05 1.46 1.88 
Dx 0.0821L DY 0.0493 
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Each point on the scatter diagram is a two-digit speci- 
men number, with the lower left-hand corner of each two- 


digit cell representing the actual point. These specimen ° 


numbers are not the specimen numbers read in with the 
data, but rather are numbers, from 1 to n, which are as- 
signed to n specimens in the data deck in the order in which 
they appear in the deck. The arbitrary set of minimum 
values is also assigned a specimen number, which is n + 1, 
and so also is the set of maximum values, the specimen 
number of which is n + 2. Where two or more points co- 
incide, a double letter is printed, with the specimen numbers 
of the coincident points printed after the same double 
letter on the right side of the diagram. 

In the sample output, specimens No, 13 and No. 14 repre- 
sent, respectively, the plots of the specified minimum and 
maximum for the plotted variables. 

The end points of the reduced major axis for the first 
diagram (DF versus AM) are plotted as follows: The 
y-intercept (the y-coordinate of the point on the reduced 
major axis where x=0), given as BASE 1! in the output 
of DASAN, is — 1.884. Because this is in the scale units of 
the diagram, it is only necessary to count down —1.9 dots 
along an imaginary vertical line through the lower left 
corner of the two-digit cell occupied by specimen No. 13, 
which represents the origin. Similarly, the point at which 
the reduced major axis crosses the right-hand margin of 
the diagram is indicated in the DASAN output by TOP 1, 
which is 2.659. This is plotted by counting 2.7 dots down- 
ward along a vertical line through the lower left corner 
of the two-digit cell occupied by specimen No. 14, the 
maximal point. 

The values DX and DY, printed in the lower right corner 
of each output page, are, respectively, the values between 
each dot on the y-axis and between every second dot on 


the x-axis. It will be seen that there are 49 such scale-units 
on the y-axis and 35 on the x-axis between the minimal 
and maximal points. 

The regression of y on x is plotted in the same manner, 
but the end points of the regression of x on y are given in 
terms of units on the x-axis rather than on the y-axis. 
Thus, TOP 3 in the output of DASAN (on the last page 
of the output listed below) refers to the number of x-axis 
units to the left (—) or right (+) of point No. 14, and 
BASE 3 refers to the number of x-axis units to the left (— ) 
or right (+) of point No. 13. 

The second diagram shown is a double log plot (LOG 
LO versus LOG AM) in which the coordinates of the speci- 
fied minimal point are (—1.0, —1.0). Here the basal end 
points of the reduced major axis and the regression of 
y on x have been measured along a vertical line through 
the dot (drawn for the purpose of illustration) which rep- 
resents the coordinates (0,0). The basal end point of the 
regression of x on y is measured along a horizontal line 
through the same dot with measurement in units of the 
x-axis. The upper end points of all three lines are plotted 
as in the preceding example. 

It will be found that artificially constructed scales match- 
ing the spacing of units on the axes will increase the speed 
and accuracy of line-plotting. When only a small number 
of points are plotted, there may be an apparent discordance 
between the plotted points (specimens) and the hand-drawn 
line of best fit. This is because the two-digit plots are only 
approximations, the accuracy of which is determined by 
the print-spacing of the printer, whereas the regression lines 
and reduced major axis, which are plotted by hand, are 
limited only by the accuracy of the hand-plotting operation 
and by the data themselves. 


Error Messages in VPLOT 


Only one error message has been written into VPLOT, 
which reads NEGATIVE VALUE IN INPUT. It is printed 
out, and the run is terminated, whenever a negative value 


appears where not permitted on a control card or in the 
data where this value is to be transformed to a logarithm. 


Availability of Program Decks 


Duplicates of the FORTRAN II source decks may be ob- 
tained through the Smithsonian Institution by writing to 


the author care of the Department of Paleobiology, Smith- 
sonian Institution, Washington, D.C. 20560. 
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